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ABSTRACT 


A nondivergent, barotropic analytical model to predict steady tropical cyclone (TC) 
propagation relative to the large-scale environment is developed in terms of a “seif- 
advection” process in which the TC is advected by an azimuthal wavenumber one gyre 
flow that results from TC-environment interaction. The model is comprehensive in that 
it includes the first-order effects of all of the dynamical influences that are presently un- 
derstood to be important to barotropic propagation: gradients of planetary and envi- 
ronmental vorticity, changes in TC wind structure, and environmental windshear. An 
unforced version of the model is used to show that angular windshear in the symmetric 
TC circulation acts to damp perturbations from axisvmmetry by tilting the perturbations 
downshear. The resultant transfer of kinetic energy from perturbation to symmetric 
circulation thus tends to restore axisymmetry. Thus, steady propagation of TC-like 
barotropic vortices is a manifestation of a stable response to asymmetric forcing. To 
predict both the asymmetric gyre flow and the propagation it induces, the forced 
Barotropic Self-Advection Model (BS.AM) is closed by seeking a particular pattern in the 
Vorticity tendencies of the TC-environmental interaction flow. For realistic combina- 
tions of environmental vorticity gradients and linear windshear, the BSAM predicts 
propagation speeds and directions that are consistent with TC propagation character- 
istics obseryed in composite data. The capability of the BSAM to account for variable 
TC structure is used to show that errors in determining TC outer wind strength of +1 
m s can result in an $5 km forecast error at 48 h. Finally, and most importantly, the 
capability of the BSAM to initialize a barotropic numerical model so that quasi-steady 
TC propagation occurs almost immediately is demonstrated for several simple dynamical 


situations. 
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Response of the model-predicted propagation velocity (columns 5,6) 
to four combinations of the parameters 6 (column 2) and y 
(colunmrs Mas defined by (SMO andmaeG) respectively. .0............ 86 


TC propagation velocity (columins 6,7) predicted by the theoretical 

model! and the numerical model (in parentheses) of Chan and 

Wilhams (1987) for three values of maximum symmetric wind 

(column 2). The analytic and piecewise-analvtic curves used by 

the numerical] and theoretical models respectively are shown in 
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Theoretically-predicted propagation velocity (columns 3.4) for 
four values of uniform environmental vorticity (column 2)............. nS) 


Theoretically-predicted propagation velocities (columns 5,6) in 

response to environmental parameter combinations (columns 2-4) 

representing TC locations south (Case 1) and north (Case 2) 

of the subtropical ridge during the western North Pacific 
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TC propagation velocities (columns 5,6) predicted by the BSAM 

for piecewise-analvtic wind profiles (see Figs. 5.2a-5.4a) that 

underestimate, approximate and overestimate the outer wind 

strength of the composite pressure-averaged typhoon (Cases 1-3 

respectively). In cach case, a quiescent environment has been 

assumed. The wind speed at the transition radius R; (column 3) 
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Original and converted composite TC motion data for the western 
North Pacific region. Column heading meanings: Vp, 1s the speed 

of the steering flow component parallel to the direction of the 
composite TC minus the speed of the TC: DD is the difference 
between the direction of TC motion and the steering flow: V. and D, 
are the speed and direction of motion of the TC respectively: 

and V, and D, are the speed and direction of the steering flow 
respectively. The data in columns V,,, and DD are taken directly 
from Chan and Gray (1982), and the data in columns V,; and D, are 
taken directly from George and Gray (1976), The data in columns V, 
and D, have been computed as described in the text. Directions 

are measured clockwise from North and the data in the last four 
columins are relative to a reference frame fixed to the surface 
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Analogous to Table A.1 for the Australian-Southwest Pacific 

region. The column headings DD, V., Dc, V, and D, have the 

same meanings as in Table A.1] and SD is the speed difference 

between the composite TC and steering. The data in columns SD, 

DD, V, and D, are taken directly from Holland (1984), except 

that the steering flow directions are measured clockwise from 

North, The data in columns V, and D; have been computed as 
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Vector differences of composite TC motion minus steering for 

the (a) latitude. (b) direction. (c) speed and (d) intensity 
stratifications of Chan and Gray (1982), and the (e) direction and 
(f) recurvature stratifications of Holland (1984). The vector 
identification labels correspond to those in column | of Tables A.1 


Relationship between the position vectors of an arbitrary point 

P in a reference frame fixed to the surface of the earth (unprimed 

variables), and in a reference frame moving at velocity C(t) 
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Hypothetical streamfunction patterns (solid, positive; dashed, 

negative) for flows that are purely azimuthal wavenumber: (a) one; 

(6) two; and (c) three. Arrows indicate general direction of the flow 
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Schematic portrayal of perturbation damping with time due to 
(a) a meridionally sheared Couette flow, and (b) a radially 
Slicarecl amisemmmiowmic wOrtical (OW. 2... ee 28 


Radial dependence of initial perturbation vorucity ( Z,. ¢* ) 

for convection-induced (k= 2, dot: k= 1, dash). motion-induced 
(chaindot), shear-induced (chaindash) and f-induced (solid) 
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(a)-(c) Perturbation streamfunction (solid, positive: dashed, 

negauve) for a k= 1 convection-induced asymmetry at t=9, 1 

and 2 hours respectively. Contour interval 1s 9.6 x 10? nr/s. 

(d)-(f) Same as (a)-(c) except for a k= 2 convection-induced asym- 
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(a) Decay of perturbation kinetic energy (normalized bv initial 

value) with time for convection-induced (k= 2, dot; k= 1. dash), 
miotion-induced (chaindot) and shear-induced (chaindash) asvmme- 

tries. (b) Same as (a). except for a f-induced (solid) and a modified 
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(a)-(c) Perturbation plus symmetric streamfunction for a f-induced 

asymmetry at t=0, $ and 36 hours respectively. Contour interval 

is 7.4 x 10° »7/s. (d)-(f} Same as (a)-(c) except showing just the 

perturbation streamfunction (solid. positive: dashed, negative), 
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44 


4.5 
4.6 


Change in symmetric vortex windspeed with time (see legend) 

as a function of radius due to a convergence of momentum flux 

associated with a 8-induced asymmetry. The windspeed of the 

mean symmetric vortex is 50 ms at the radius of maximum winds 

and decreases with inverse dependence on radius to 5 m/s at ten 

tinies the radius of miaximunt winds. ... 20). Bee ns Sas aes 


(a)-(c) Perturbation vorticity (solid, cyclonic; dashed, anticyclonic) 
for a B-induced asvmmetrv at t=0, 1 and 2 hours respectively. 
Contour inter.alis Gabe IVs Ar ee eee cn ee, 48 


An illustration of the model domain and three annular subdomains. 

The inner and outer boundaries of the model domain are denoted 

by R, and R; respectively. The interface between the inner (n= 1) 

and outer (n= 2) annulus of the Self-advection Region is denoted 

by R,. The interface between the Dispersion Region annulus (n= 3) 

and the inner annulae is denoted by R, and corresponds to the 
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(a) Radial profiles of the analytic TC windfield for the Fiorino 

and Elsberrv (1989) basic vortex (dashed) and the piecewise-analytic 

function used (solid) to approximate it. Parameters that define 

the piecewise-analvtic function according to (4.6) and (4.16) are 

shown in the inset. Vertical dotted lines (left to right) correspond 

to the radial boundaries interfaces R,. R, and R,(=R,). The 

chain-dashed curve represents the function f,,.7?. (b) The analytic 

and piecewise-analvtic vorticity gradient profiles associated with 

the Willd prCilesan (alee, re ee en rene a Poe em: ee | 


Model-predicted wavenumber one gyre streamfunction (solid, 
positive; dashed, negative) using the piecewise-defined wind profile 
of Fig. 4.2a, and the parameter specifications: Ty.=40 m/s, Rye = 


100 km, f.= 2x 10-8 rts, 6 = 90, Rp = 1.5, R,= 100, C.= 2.65 ms 

and »=120, The contour interval is 2 x 10°172/s, and only the inner 

2400x2400 of the domain is shown to correspond to the illus- 
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Numerically predicted asymmetric streamfunction at 48 h due 
to f-induced distortion of a TC wind profile initially defined as 
in Fig. 4.2a. Contour interval is 2 x 10'n7?/s, and the distance 


between axis tick marks is 40 km (Fiorine and Elsberry 1989),..0....... 67 
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(a)-(b) Vector diagrams showing direction and magnitude of £, 

tendency maxima that have been scaled by — 6{,/ér at (a) r= 2 and 

(b) r= respectively. The vectors represent the terms of (4.24) 

(see key above), and were computed using the piecewise-defined 

symmetric TC wind profile of Fig. 4.2a and the wavenumber one 

evre solution of Fig. 4.5. (c) Radial profiles of the amplitude (solid) 

and phase (dashed) of {, corresponding to the streamfunction 

field of Fig. 4.5. The amplitude curve has been scaled by the 

amplitude of ¢, at r=0.2. Vertical dotted lines (left to right) 
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(a)-(c) As in Fig. 4.7, except using the analytical model solution 
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Vector differences between the TC propagation and an averaged 
interaction flow velocity (0-300 km) at various times (see legend) 
during a 120 h integration of a barotropic numerical model (Fiorino 
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As in Fig. 4.7, except using the analytical mode] solution 
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Radial profile of streamfunction amplitude (solid) and phase 

(dashed) corresponding to vorticity profile shown in Fig. 4.10c. 

The amplitude curve has been normalized by the maximum 
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As in Fig. 4.5, except using the parameter specifications of Test 4 
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(a)-(b) As in Fig. 4.2, except that the parameters (inset) if the 
piecewise-analvtic TC wind structure have been changed so that 

the associated vorticity gradient is no longer continuous and 

egualtouhe analvtic value at the transition radius (y=4.75)............. 89 


As in Fig. 4.2. except using Chan and Williams (1987) analvtic 


WE ne weeoiile paramicters of 1 ).—20 ms, K,.— 100 km and b=1.0...... o] 
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4.24 


5.4 


Theoretically-predicted wavenumber one streamfunction fields 
for (a) Case 1 and (b) Case 4 of Table 4.3. Contour interval 
18 2% OES. ee ees sn eae ee CRE Gee ae a7 


Idealized planetary and environmental vorticity gradients for 

TC positions in an anticyclonic vorticity region (a) south and 

(b) north of the subtropical ridge in the lower to middle troposphere 

of the western North Pacific TC region. The units of the vorticity 

gradients are 10-"(7is)-", ame the units of 4, ane 10ers 98 


As in Fig. 4.18, except that the piecewise-defined wind parameters 
(see inset) have been recalculated to account for B,~= 3.5x 10“ ms"! 
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(a)-(b) Vorticity tendency vector diagrams, and (c) vorticity 
amplitude and phase profiles as in Fig. 4.7, except for the model 
solutiomHorn Gases) Paper dian cn se le eae a as 4 101 
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Model-predicted wavenumber one streamfunction fields as in 
Fig. 4.3, except for (a) Case 1 and (b) Case 2 of Table 4.4. 
Contourintenr val is 92a Gia) sree a ern rss 4 soc ane 103 


(a) Radial profiles of Fiorino and Elsberrv’s “basic vortex” 

(dashed), the composite surface tangential winds of large Atlantic 

hurricanes (solid squares; Merrill 1984) and the tangential winds 

of the composite pressure-averaged tvphoon (open squares; Frank 

1977). Here T',.=35 m’'s and R,.= 100 km. (b) Analytic approx- 

imations to the radial profile of the composite pressure-averaged 

tvphoon (open squares) as defined bv (4.22) with b= 0.6 (dashed) 

and (3,1) (solid: paraineterssimimset jas ieren| yom OMG eee 5 107 


(a) Radial wind profiles of the composite pressure-averaged typhoon 

(dashed) as approximated by (5.1) in Fig. 5.1b and a piecewise- 

analytic profile (parameters in inset) that overestimates outer 

wind strength. (6) Piecewise-analvtic and analvtic vorticity gra- 

dient profiles corresponding to the winds in (a). In both cases, Vy.= 
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As in Fig. 5.2, except for a piecewise-analytic wind profile that 

closely approximates the outer wind strength of the composite 
pressure-averaged tvphoon in the 300 to $00 km radius “critical 
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As in Fig. 5.2, except for a piecewise-analvtic wind profile that 
underestimates the outer wind strength of the composite pressure- 
averaged typhOOlii cs jcc sceeccesee te ox ae eco ays on 11] 
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Wavenumber one gvre streamfunction (dashed, negative) associated 
with Table 5.1 Case 2. Contour interval is 2.0 x 10°72/s. Dotted 
circles have radii of 2°, 4°, 6°, $° and 10° lat. The middle circle 
corresponds to the center of the 5°-7° lat. radius annulus typically 
used to compute steering from composite data (e.g., Chan and 


Gray 1982, Holland 1984). Theparameter.&..= 100'km. ........... 


(a) BSAM-predicted propagation (C) for the composite pressure- 
averaged typhoon and contributions to environmental steering bv 
the associated BSAM-generated wavenumber one gyre flow eval- 
uated around the 2.4.6.8 and 10° lat. radius circles of Fig. 5.5. 

(b) Differences between composite pressure-averaged typhoon 
propagation (C) and the wavenumber one gvre contributions to 
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As in Fig. 3.5, except for including in steering flow computation 
only the winds in the 90° arcs to the right and left of the plane- 
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As in Fig. 5.6, except for vectors calculated based on the arcs 
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f-induced propagation tracks and 6-hour positions predicted bv 
the numerical model of Chan and Williams (1987) for an initial 
TC wind profile corresponding to Table 4.2 Case 2 from 0-48 h 
(solid circles), and for the period 36-84 h (crossed circles) except 
translated so that the 36 h position corresponds to the initial 
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Wavenumber one gvre streamfunction patterns generated by the 
BSAM for the ssmmetric 1C of Table 4.2 Case 2 for (a) Chan 
and Williams propagation velocity of 2.8 ms at 330° (x= 120°), 
and (b) a BSAM-predicted propagation velocity of 2.65 m § at 

%= 132°. Areal extent of the figures corresponds to the domain 

of the numerical model. and the streamfunctions have been linearly 


adjusted to zero within 20 gridpoints of the domain boundaries....... 


The f#-induced propagation tracks and 6-hour positions predicted 
by the numerical model initialized with the symmetric TC as in 
Fig. 5.9, except including the gvre structures of Fig. 5.10a (solid 
circles) and Fig. 3.10b (open circles). For comparison, the trans- 


lated 36-84 h track (crossed circles) from Fig. 5.9 is also shown. ...... 


Tracks and 6-hour positions predicted bv the numerical model for 
f-induced TC propagation in an zonal current with linear anti- 
cualpnic shear of Z.. = —0ix4 10s using the svmmetric 

initial TC wind profile of Fig. 5.9 with (open circles) and without 


(solidicirclesy) BSAM-eenerated wavenumber One gires............. 


i iontie@ eae except ior cyclonic shear of Z7,. = 3.0x ls!) ...... 
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I. INTRODUCTION 


It is Well known that tropical cyclone (TC) motion persistently differs from the mean 
advective effect (i.e., the “steering”) of the large-scale environment (e.g., George and 
Gray 1976; Chan and Grav 1982). The precise speed and direction of this velocity dif- 
ference (hereafter referred to as “propagation”) depends to some extent on the particular 
steering flow definition employed, but is generally westward and poleward in direction 
and I to 2.5 ms in magmitude (Carr and Elsberry 1989). Although this propagation 
velocity 1s typically smaller than the advection by environmental steering. it can induce 
sigmficant TC track errors (2 ms = 170 km d) if not properly accounted for in dyvnam- 
ical models. Interestingly. it is well known that numerical model forecasts of TC tracks 
generally have less skill than CLIPER (CLImatology and PERsistance) over the first 
24-36 h, and that the error is primarily an underestimation of TC speed (e.g., Neumann 
and Pelissicr 1981: their Figs. 4.8 respectively). Since TC’s track generally westward and 
poleward throughout much of their lifetimes, such statistics provide circumstantial evi- 
dence that dynamical forecast models do not as vet properly account for TC propagation 
unul about 1 day into the numerical integration. It will be shown that the initial 
slowness of numerical models mav result from initializing the model with a symmetric 
TC structure ther lacks any propagation-indJucing asymmetries caused bv 
TC-entironment interaction. Although past modeling experience indicates that both 
barotropic and baroclinic processes can cause TC propagation asymmetries, the general 
goal of this research is to develop an mmproved theoretical model for understanding and 


predicting the barotropic contributions to TC propagation asymmetries. 


A. MODELING REVIEW 
1. Status of basic understanding 
a. f-induced propagation 

The propagation of an iniually symmetric, nondivergent barotropic (NDBT) 
Vortex in a quiescent environment due to the influence of planetary rotation, as ap- 
proximated by a £-plane. provides the most basic dynamical scenario for beginning a 
review of TC motion theorv. The northwestward propagation at 2-3 ms noted by Chan 
and William (1987; hereafter CW) in a recent numerical study of this problem is gener- 
allv consistent with the propagation observed in previous numerical studies (Anthes and 


Hoke 1973: Kitade 1981; DeMaria 1983, hereafter DM). The authors of the previous 


numerical studies attributed their results to the traditional linear theories of barotropic 
vortex propagation advocated in various forms (e¢.g., Rossby 1939, 1948: Kasahara 1957; 
Kasahara and Platzmann 1963; Holland 1983). An interesting aspect of these theoretical 
predictions is that the westward and poleward components of the propagation are asso- 
ciated with different linear mechanisms, e.g., Rossby wave propagation by the TC to the 
west as opposed to poleward acceleration due to net Coriolis force on the TC. In con- 
trast, CW concluded that only the nonlinear “self-advection” (i.e., the TC advecting its 
own Vorticity) process associated with the development ofa horizontal asymmetry in the 
TC windfield was responsible for both the westward and poleward components of the 
propagation. 

The basis for this important distinction can be demonstrated by partitioning 
the vortex flow into a steady symmetric (S) component and an asymmetric (A) compo- 
nent that is generally small by comparison. In a coordinate system centered at any in- 
stant on the symmetric component, the propagation process in a quiescent environment 


is Approximately described by 





et = = V ye Vig — V5-Vey - V5- Wf, (1.1) 
(a) (b) (©) (d) 


where term (a) represents the vorticity tendency associated with motion and or dis- 
tortion of the TC due to the combined effect of the two dominant self-advection terms 
(b and c) and the advection of planetary vorticity (d). Linear theories regard the west- 
ward component of propagation as a Rossby wave-like propagation associated with the 
east-west dipole of vorticity tendency generated by term (d), and either eliminate self- 
advection by requiring the TC to be symmetric at all times (e.g., Kasahara and Platzman 
1963), or assume that self-advection is implicit in computations of observed steering 
flows (e.g., Ilolland 1983, 1984). When the nonlinear self-advection term was omitted 
from their model. CW noted that advection of planetary vorticity by the TC 1s a strong 
dispersive. distorting influence, but by itself generates negligible motion of the TC center 
(=0.3 ms). Thus, they concluded that the linear process generates the asymmetric flow 
essential to the existence of the self-advection process, but the self-advective interaction 
of the symmetric and asymmetric TC components actually causes the propagation. 
Tiormo (1987) and Fiormo and Elsberry (1989; hereafter FE) clarified the 


individual roles of the two primary self-advection processes (terms b and c) via a diag- 


to 


nostic technique that extracts V, and V, from a numerical model similar to that used 
by CW. Their analvsis showed that the structure of V, is predominantly an azimuthal 
wavenumber one asymmetry that consists of two counterrotating gvres (hereafter re- 
ferred to as “Wavenumber one gvres”) that produce a near-uniform flow across the cen- 
tral region of the TC. The average velocity of the uniform flow is very nearly equal to 
the TC propagation vector, which indicates that NDBT -induced propagation is es- 
sentially an advection of symmetric TC vorticity by the asymmetric flow component of 
the IC, i.e., term (b) of (1.1). This is somewhat surprising since simple scale arguments 
suggest that the advection of asymmetric vorticity by the symmetric tangential wind 
(term c) should be equally important to the propagation process. However, FE demon- 
strated that streamfunction tendencies associated with term (c) nearly cancel those as- 
sociated with term (d) after a quasi-steady state is achieved. Thus, term (c) acts 
primarily to linut the growth of the wavenumber one gyres and the resultant propagation 
to nearly steady values after an initial adjustment period of approximately 24-48 h. 

A theoretical model that accurately predicts f-induced gvre structure and 
the associated propagation based on this recently identified self-advection process has 
vet to be developed. and this is the principal objective of this research. The successful 
development of such a model that confirms the numerical work of CW and FE will be 
an important step in establishing the nonlinear self-advection principle as the mechanism 
responsible for barotropic TC propagation as opposed to the traditional linear Rossby 
wave arguments. 

‘Since advection of planetary vorticity by a vortical flow represents essen- 
tially steady asymmetric forcing, FE’s results suggest that advection by the svmmetric 
tangential wind tends to stabilize a barotropic approximation of a TC to such forcing. 
Nonhnear advection is necessary to permit dispersion-resistant movement of model 
Vortices in general (e.g.. \icWilhams and Flierl 1979; Mied and Lindemann 1979), or 
nondispersive movement in the case of very specialized entities known as either modons 
or solitary eddv solutions (see Flierl er al. 1980). However, the precise mechanism by 
Which nonlinear advection stabilizes a vortex of arbitrary radial structure has not been 
identified. Because steady barotropic vortex propagation is essenually a manifestation 
of this stability to asymmetric forcing. a thorough analvsis of this process is important 
to improved understanding of barotropic TC propagation, and thus is another objective 
of this rescarch. 

It is important to briefly address the question whether barotropic vortex 


propagation stability is a sufficiently important process compared to potential baroclinic 


influences to warrant the present barotropic modeling effort. Although Holland (1983) 
has suggested that inertial stability accounts for near axisymmetry in the inner 300 km 
of a TC, inertial stabilization cannot occur in NDBT or quasi-geostrophic model. The 
qualitatitive similarity of TC propagation tracks in simple baroclinic models (e.g., 
Madala and Piacsek 1975; Kitade 1980) to barotropic results provides at least 
circumstantial evidence that baroclinic TC stability’/propagation is more a modification 
to, as opposed to being fundamentally different from, the barotropic phenomenon. 
Thus, it is reasonable to assume that TC propagation/stabilitv is fundamentally a 
vorticity advection process that should be fully understood within a barotropic context 
before additional baroclinic complexities are added. 
b. Influence of TC structure 

By using an initial vortex specification that permits independent variations 
of the inner and outer tangential wind structure, FE clearly showed TC propagation 1s 
strongly dependent on outer wind changes (strength), but is virtually independent of in- 
ner wind changes (intensity) as predicted by earlier researchers (Holland 1983; DM). In 
particular, an increase of only a few m s in the initial tangential wind between 300 and 
800 km results in a faster and more westward TC propagation that can significantly alter 
the f-induced TC track. Although the initial value of TC total relative angular mo- 
mentum (RAM) has been cited as having a potentially important influence on TC tracks 
(DeMaria 1987), Shapiro and Oovama (1989) have recently shown that f-induced mo- 
tion does not depend on initial TC RAM in any well-defined way. Since RAM 1s es- 
sentially an integral (and therefore not unique) measure of TC strength. TC propagation 
variations due to changes in RAM are likely a manifestation of the dependence on TC 
strength described by FE. 

In view of the above, it is clear that properly initializing dynamical forecast 
models for TC strength, and ensuring that the model dynamucs’physics maintain an ac- 
curate TC outer wind profile throughout the integration are important areas for ongoing 
research. Without aircraft, improved measurements of TC structure for initializing op- 
erational forecast models can only be achieved by advancement in remote sensing ca- 
pabilitics. Designing utilizing remote sensors to better detect TC strength and choosing 
mode] parameterizations to best maintain a desired TC structure both depend on im- 
proved understanding of the nature and sensitivity of TC propagation dependence on 
outer wind strength. Thus, an additional objective of this research is to identify the dv- 


namical basis for this dependence and provide at least a preliminary assessment of the 


Wind measurement accuracy required to adequately account for the influence of TC 
outer wind strength on the barotropic propagation process. 
c. Influence of divergence - 

Anthes and Hoke (1975) noted significant differences in the B-induced TC 
motion tracks predicted by nondivergent and divergent (5 km depth) barotropic numer- 
ical models, which they attributed to divergence-induced slowing of Rossby wave phase 
speeds. Such a result would suggest that nondivergent forecast models such as 
SANBAR (Sanders es al. 1975) fail to account for a significant propagation process. 
In a model comparison similar to Anthes and Hoke, Shapiro and Ooyama (1989) found 
nearly identical propagation tracks, even though substantially more divergence was in- 
cluded (1 km depth). Thus, Shapiro and Oovama concluded that the earlier result was 
due to dissinular initial TC structures that were emploved in the divergent and nondi- 
vergent models. This result is important for clarifving the minimal role divergence plays 
in barotropic TC propagation, and also supports the choice of a NDBT dynamical 
framework for the theoretical model emploved in this study. In addition, their result re- 
emphasizes the important influence of initial TC specification on the accuracy of TC 
tracks predicted by a barotropic model. 

d. Environmentally-induced propagation 

Barotropic models also have been used to study the influence of horizontally 

nonuniform environmental winds on TC propagation relative to steering. The additional 


influences on the propagation process may be described by rewriting (1.1) as 
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in which the environmental (E) processes V,+V(f+¢,) and €¢,/ér have been omitted 
since thev relate to the evolution of the steering flow. Terms (a). (b) and (c) of (1.2) are 
the same as in (1.1). However, term (d) of (1.2) now includes advection of environmental 
vorticity by the TC, and a new term (e) accounts for both steering and shearing of the 
symmetric TC and wavenumber one gvres by the environmental wind. 

With the addition of variable environmental winds, it becomes essentially a 
matter of convention whether V, should be associated with the symmetric TC (V) in 
which case the term “self-advection” is still meaningful. or whether V, should be re- 


garded as a modification to V,;. and thus a contribution to steering. The comparison of 


observational data with theory presented in Section C will provide strong evidence to 
suggest that V, will be manifested as propagation rather than a contribution to con- 
ventionally computed environmental steering. Thus, V, will continue to be referred to 
as a self-advection flow associated with TC-environment interaction, and will be re- 
garded as distinct from the steering associated with V;. 

In a NDBT numerical study with steady zonal winds that varied only with 
latitude, DM noted that the poleward component of propagation associated with a TC 
in a poleward (equatorward) gradient of environmental vorticity increased (decreased) 
the northward propagation component relative to f alone. DM attributed his results to 
the linear theory of Kasahara and Platzman (1963), which predicts that horizontal vanri- 
ability of the environmental winds will induce TC propagation in the direction of and 
90° to the left of the large-scale gradient of absolute vorticity. This theorv is essentially 
a direct extension of Rossby’s westward propagation (1939) and poleward acceleration 
(1948) of TC’s due to f alone. The presence of £, in term (d) of (1.2) would make 
Kasahara and Platzman’s hypothesis intuitively plausible. However, the recent demon- 
stration by CW that linear advection of planetary vorticity only distorts the TC must 
also apply to the advection of environmental vorticity by the TC. DM’s results clearly 
indicate that environmental vorticity gradients alter barotropic TC propagation. This 
suggests that Kasahara and Platzmann’s theorv of a linear absolute vorticity-induced 
propagation due to the combined effects of # and ¢, can be recast within the context of 
the nonlinear self-advection mechanism recently identified. Thus, another objective of 
this research will be to determine the extent to which a self-advection model for 
f-induced propagation can also accomodate the presence of a horizontally variable 
windfield. 

Anv analogy between the propagation-inducing effects of f and {, is com- 
plicated by the additional distorting influence of environmental wind shear in term (e) 
of (1.2). In a barotropic model that includes # and an environmental current with only 
linear shear, Chan and Williams (1989) found that cyclonic (anticyclonic) shear induces 
a westward (eastward) curving TC propagation relative to the tracks induced by f alone. 
Thev also noted that the orientation of the wavenumber one gvres rotates consistently 
with the propagation changes. which suggests that this process is due to shearing of gyre 
vorticity by the rotational part of the environmental current. Although curving TC 
tracks are also evident in the results of DM who used a sinusoidal windfield, the effect 


is much smaller than that seen by Chan and Wiliams. Because of the potentially im- 


portant impact of environmental shear on TC propagation direction, modeling of this 
linear shear process Will be included in this research. 
2. Status of propagation/gyre prediction models = 
a. Recent results 

The preceding discussion of the self-advection mechanism for barotropic 
TC propagation has been based on diagnostic analyses of numerical modeling results. 
Theoretical models of TC motion that explicitly predict wave number one gyre structure 
and associated propagation based on the self-advection concept are also being devel- 
oped. The limited success of such models as discussed below has provided motivation 
for this research. 

Given the intense advective influence of the TC’s tangential winds, the 
quasi-steady nature of the £-induced wavenumber one gyre observed by FE suggests that 
It is a standing azimuthal wave. Such a wave would propagate clockwise in a Rossby 
wave-like manner on the basic state vorticity gradient represented by the symmetric 
component of the TC via term (b) of (1.1). On the basis of such a hypothesis and an 
assumption as in (1.1) and (1.2) that V,<< V, over a large area, Willoughby (1988) 
formulated a linear shallow model that explicitly predicts the structure of such “vorticity 
Waves” as well as the propagation they induce. Ile solved the potential vorticity and 
divergence equations in a polar coordinate system moving with the TC for a range of 
possible propagation velocities, and assumed that the correct propagation velocity would 
munimize the Lagrangian of the system as required by Hainilton’s Principle from vanri- 
ational calculus. By incorporating a rotating mass source-sink pair to approximate 
asymmetric evewall convection. he obtained a vorticity wave angular frequency that 
agrees well with that of small-scale trochoidal oscillations observed in TC tracks. This 
result suggests that TC track oscillations are also explainable in terms of a barotropic 
self-advection process associated with asymmetric convection, as opposed to the tradi- 
tional Magnus force theorv (Kuo 1950; Yeh 1950). 

Willoughby’s model failed to predict a correct propagation velocity on a 
f-plane (quiescent environment), presumably due to a barotropic instability studied by 
Peng and Williams (1989) in a nondivergent version of the model. Peng and Williams 
showed that existence of the instability depends on the vorticity gradient sign reversal 
of the TC (necessary condition) and sufficient model resolution to resolve the small scale 
of the instability (radius of maximum winds). Thus, they suggest that this linear insta- 
bility is dynamically based, but is highly damped in full numerical models due to a 


combination of nonlinear interaction and coarse horizontal resolutionn. Jlowever, ana- 


=) 


lytical gyre propagation models such as the one discussed in the next paragraph and the 
one developed herein do not experience the instability. Further analysis of this process 
is necessary to resolve the apparent paradox, but this issue is not considered within the 
scope of the present research. 

Smith ef al. (1989) have devised an approximate analytical model for the 
B-induced gvre structure by ignoring vortex propagation and computing the 
streamfunction associated with a redistribution of the initial absolute vorticity by the 
symmetric TC. The inner gyre uniform flow predicted by the model closely approxi- 
mates the propagation speeds predicted by nonlinear numerical models up to about 24 
h. Smith (personal communication, 1989) is presently modifying the model to extend the 
period of accurate propagation prediction by empirically accounting for the movement 
ofthe 1 

b. Application 

A natural application for a barotropic TC propagation model would be as 
part of a track prediction scheme that computes total TC motion at any time by 
adding: i) a steering velocity obtained from a numerical model for the evolution of the 
large-scale environmental winds; and ii) a propagation velocity derived from the theo- 
retical model based on the gradient of absolute (or potential) vorticity of the numer- 
ically generated environmental winds. Early efforts based on linear propagation theory 
(Kasahara 1957; Kasahara and Platzman 1963) included the cross-vorticitv gradient 
propagation effect, but ignored the hypothesized up-vorticity gradient acceleration of the 
TC. In fact, this supposed acceleration due to a net force on “the vortex” has never been 
transformed into a calculated propagation velocity. A recent track prediction model 
developed by Holland (1983: 1984) computes westward propagation due to # similar to 
the earlier efforts, but also includes a northward propagation effect by empiricallv ac- 
counting for net convergence into the TC. That these models have shown only limited 
forecast skill is not surprising in light of the recent self-advection theory of propagation, 
and the skill that they do possess is probably due to some portion of the actual nonlinear 
self-advection phenomenon being parameterized by the linear models. 

The present unavailability of an accurate TC propagation model based on 
nonlinear self-advection has precluded any implementation into an “advection + prop- 
agation’ track prediction model analogous to the linear models. Although the present 
model may prove suitable for this type of endeavor, development and testing of such a 
model is outside the scope of this research. Unlike the linear propagation models that 


are insensitive to small deviations of the 1C from axisymmetry, self-advection models 


depend fundamentally on such asymmetries, and more importantly, actually predict their 
structure (e.¢., Willoughby 1988; Smith ez a/. 1989). Accurate predictions of such gyre 
structures would be ideally suited for initializing barotropic numerical models such as 
SANBAR (Sanders er al. 1975) so that quasi-steady propagation occurs immediately at 
the start of the integration rather than after the 24-48 h typically observed for an initially 
symmetric TC (e.g., DM: CW; FE). This dissertation will conclude with a proof-of- 
concept demonstration of how the gyre propagation model developed here can be used 


to improve numerical model initialization. 


B. OBSERVATIONAL EVIDENCE FOR PROPAGATION 
Given the manv barotropic models of TC propagation in the literature, it is natural 
to ask whether observational evidence exists to support the predictions of the theoretical 
or numerical models. Unfortunately, data deficiencies have made computation of prop- 
agation for an individual TC difficult. However, studies that average the data from 
many TC’s to produce a “composite” TC are available for this purpose. The propa- 
gation vectors in Fig. 1.1 from Carr and Elsberry (1989) are based on data from two such 
studies (George and Gray 1976; Chan and Gray 1982). A description of the process by 
which the vectors were obtained 1s given in Appendix A. 
}. Comparison with modeling results 
The propagation vectors in Fig. 1.1 exhibit a number of interesting properties 
that strongly resemble the TC propagation in the numerical models cited above. Except 
for the anomalous “after recurvature” vector (Fig. 1.1f), the vectors have magnitudes 
ranging from 1.0 to 2.5 ms and directions that tend to be westward and poleward in 
both hemispheres, which is consistent with the numerical and analytical results previ- 
ouslx cited. In addition, the rotation of the propagation vector direction from west- 
southwestward for westward moving 1C’s to northwestward for eastward moving TC’s 
in the direction stratification (Fig. 1.1b) is consistent with DM’s numerical results. DM 
showed that the change in the direction of the environmental vorticity gradient from 
poleward on the poleward side of the subtropical ridge to equatorward on the 
equatorward side caused a decrease in the meridional component of TC propagation 
similar to that in Fig. 1.1b. Finally. the propagation vectors in the imtensity 
stratification (Fig. 1.1d) have a direct dependence on TC intensity. As discussed earlier, 
the modeling studies of DM and FI: demonstrate that a nondivergent barotropic pre- 
diction of TC propagation duc to f is independent of TC intensity. but is well correlated 


with outer wind strength (see Merrill 1984 for tvpical definitions of strength and inten- 





Fig. 1.1 Vector differences of composite TC motion minus steering for the (a) latitude, 
(b) direction, (c) speed and (d) intensity stratifications of Chan and Gray (1982), and the 
(e) direction and (f) recurvature stratifications of Holland (1984). The vector identifica- 
tion labels correspond to those in column I of Tables A.] and A.2. 
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sity). Since a weak correlation exists between the intensity and strength of TC's 
(Weatherford and Gray 198$), the increase in propagation vector magnitude in Fig. 1.1d 
may be a manifestation of the numerically-predicted dependence of $-induced-propa- 
gation on TC strength. 

The results in Fig. 1.1 also contain some apparent inconsistencies. Examples 
are: 1) the significantly larger meridional components of the Northern Hemisphere vec- 
tors compared to that of the Southern Hemisphere vectors; and i) the presence of 
equatorward components in some of the propagation vectors. Such properties may be 
associated with boundary layer or baroclinic processes not considered in the barotropic 
theories. However. a possible barotropic explanation might be the presence of east-west 
vorticity gradients in the TC environment that were excluded by DM. For example, a 
large-scale Westward relative vorticity gradient is present during the summer in the 
troposphere between the anticvclone over the western North Pacific and the heat low 
over southeastern Asia. Based on DM's results, a TC vortex embedded in such a 
Vorticity gradient should have a westward, and more importantly, a southward compo- 
nent of propagation. Since the meridional gradient of environmental relative vorticity 
and f# are in opposite directions south of the Northern Hemisphere subtropical ridge, the 
zonal gradient of environmental vorticity might tend to dominate, and thus explain the 
southward component of vector W in Fig. 1.1b. In contrast, the meridional environ- 
mental vorticity gradient and f are in the same direction north of the Northern Hemi- 
sphere subtropical ridge. and thus might dominate over the influence of a zonal vorticity 
gradient in the cases of vectors N and E in Jig. 1.!b. Differences in the direction of the 
large-scale absolute vorticity gradients in the Northern and Southern Hemispheres thus 
may contribute to the hemispheric variability of the data shown here. The present 
propagation model will permit testing of such a hypothesis within a modeling context 
since the influence of an east-west environmental vorticity gradient will be included. 

Statistical influences in Fig. 1.1 also must be considered, such as: 1) ambiguities 
introduced by composite stratifications that may incorporate multiple propagation- 
inducing influences; and 11) possible random or systematic errors in the composite data 
that may be significant relative to the small size of the propagation vectors being ana- 
Ivzed (e.g.. accuracy of rawinsonde wind measurements and error in locating the TC 
center position). Random and systematic errors should be reduced as sample sizes are 


increased and observational accuracies are improved respectively. 


ul 


2. Interpretation of observed propagation vectors 

The fundamentally important issue that determines the mechanism responsible 
for observed TC propagation is whether or not a TC-interaction flow (V,) will be man- 
ifested as an essentially inseparable contribution to conventionally computed 
enVironnmental steering. If V, is included within the steering, then the observed TC 
propagation vectors must be due to some mechanism other than self-advection. at least 
within the context of the barotropic theory. Ifa self-advection flow is not accounted for 
in the computed environmental steering flows, then the observed propagation vectors 
may indeed represent the self-advection process. 

Holland (1983) assumed that V, 1s included in the steering, and he has proposed 
a linear propagation model that is consistent with the near-zonal orientation of the 
Austrahan-Southwest Pacific difference vectors (Fig. l.le-f). Holland assumes that 
inertial stability constrains the inner core of the TC to move with the outer envelope, 
which is assumed to propagate westward as a Rossby wave with a phase speed appro- 
priate to an “effective radius.” In practice, the “effective radius” parameter for a partic- 
ular storm and time is chosen to give a barotropic Rossby wave propagation speed that 
equals the observed westward component of TC propagation over a preceding time in- 
terval, Although Holland proposed low-level convergence to account for small devi- 
ations from pure westward motion, the above discussion of environmental shear in 
addition to f# is an alternative explanation. 

In contrast to Holland's assumptions, FE’s explicit illustration of the 
wavenumber one structure of the self-advection flow provides strong evidence that V, 
1s, for the most part, not included in the steering. Key aspects of this structure are that 
the central uniform flow portion of the gvres that account for the propagation is con- 
fined to radii less than 300 km from the TC center, and that the gyre flow is considerably 
weaker outside that radius. Even allowing for somewhat larger gvres in actual TC’s, this 
suggests that the uniform flow region of the interaction flow would be largely unac- 
counted for in a steering flow calculated over an annulus of 5°-7° lat. radius from the 
TC center as is the case for the data used to produce Fig. 1.1. Thus, “self-advection” 
should be manifested primarily as propagation, rather than as a contribution to con- 
ventionally calculated steering flows. Thus, Holland’s linear model may actually include 
nonlinear motion-inducing processes since the selection of an effective radius is based 
on the observed difference between TC motion and steering over a preceding time in- 


terval. 


Accurate measurement of steering near the TC center will continue to be prob- 
lematic for the foreseeable future. Thus, it seems advisable as a practical matter to 
continue to compute steering at a large scale (= 1000 km) as is now done, and to regard 
the self-advection flow as a propagation-inducing process that is distinct from steering 
by the large-scale environment. In light of the subsynoptic scale region in which the 
self-advective flow is influential, such a flow partitioning may also be a good approxi- 
mation to the long-sought scheme to uniquely separate the TC from its environment. 
Thus, this research will use the three-part partitioning scheme proposed in Elsberry 
(1986) that consists of: 1) a symmetric TC; 11) a propagation-inducing TC-environment 
interaction flow; and iii) a large-scale environmental flow component that accounts for 


TC steering. 


C. DISSERTATION OVERVIEW 
In summary, the specific objectives of this research are: 


1. identification and analvsis of the specific mechanism responsible for barotropic 
Vortex stability to asymmetric forcing, and how this stabilizing influence is related 
to TC motion; 


2. development and analvsis of a theoretical model for barotropic TC propagation due 
to both planetary and environmental vorucity effects that is based on nonlinear 
self-advection, and that accurately predicts gyre structure and associated propa- 
gation velocity relative to equivalent numerical model solutions; 


ted 


identification of the dynamical basis for the dependence of barotropic TC propa- 
gation of outer wind strength, and an assessment of the wind measurements re- 
quired to account for this phenomenon: and 


4. provide a preliminary demonstration of the viabilitv of self-advection propagation 
models for improved initialization of barotropic TC track forecasting models. 
Because of the mathematical complexity of the model to be constructed, Chapter I] 
will be devoted to preliminary development and analysis. Objectives | and 2 above will 
be addressed in Chapters III and IV respectively, and the theoretical aspect of Objective 
3 will be addressed in Chapter IV as well. In Chapter V, the model will be applied to 


satisfy the second half on Objective 3 and Objective 4. Chapter VI will conclude the 


dissertation with an overview of the modeling approach and a summary of results. 


H. PRELIMINARY MODEL DEVELOPMENT 


A. REFERENCE FRAME TRANSFORMATION 
The behavior of fluid motion in a NDBT system on a f-plane is governed by the 


conservation of absolute vorticity, which may be expressed as 


“ap ts aes (2.1) 
C= how ve (2.2) 
f=Jo + By, (2.3) 


where ¢ is the local vertical component of relative vorticity, 8 is the linearized latitudinal 
derivative of the Coriolis parameter f, and f, is the average value of fin the domain. In 
terms of Eulerian partial derivatives, (2.1) may be expressed as 


ise 


+ VME +f) a0. (2.4) 





A formal partitioning of the total fluid flow that takes advantage of the near- 
axisyinmetry of TC’s as suggested by (1.1) and (1.2) requires (2.4) to be transformed to 
a reference frame moving with the TC at a translation velocity C(z). In terms of 
cartesian coordinates, the relationship between the position vectors for an arbitrary 
point P in each reference frame is illustrated in Fig. 2.1, and is mathematically defined 
by 


R(x, y) = R’(*’, 39”) + [coer ; 23) 
0 


Where the (") symbol denotes variables in the moving reference frame. Let the (’) symbol 


also define certain derivative operations in the moving reference frame 
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Fig. 2.1 Relationship between the position vectors of an arbitrary point P in a reference 
frame fixed to the surface of the earth (unprimed variables), and in a reference frame 
moving at velocity C(t) (primed variables). 


By using the chain rule and (2.5)-(2.7), it may shown that 





eee Oey (2.9) 
Cl 

which, when substituted into (2.4), results in the desired transformed expression 

‘| ~cev oe vere ys). (2.10) 


It is important to note that a convention of leaving the dependent variables 
untransformed in (2.10) has been emploved. Although relative vorticity and the gradient 
of the Coriolis parameter are invariant with respect to the coordinate transformation, 


fluid velocity is not. Thus, the velocity transformation 
Vow + Ci. (2.11) 


obtained from d dt of (2.5) has not been used. This results in explicit advection by C in 
(2.10). but has the desirable property of avoiding an implicit dependence of the fluid flow 
on C in the moving reference frame. Such a convention is necessary to justify the use 
of the homogeneous boundary conditions that will be emploved in solving this model, 


and facilitates later use of the solutions in a numerical model initialization scheme. 


B. FLOW PARTITIONING 


As outlined in Chapter I, the total fluid flow is partitioned into 
V = Vs + VE + Ny ’ (202) 


according to the following definitions: i) V,; is a known symmetric /S); TC flow com- 
ponent that is steady in the moving reference frame; ii) V; is a known environmental 
(E) flow: and iii) V, is a predominantly asymmetric (4) flow component that repre- 
sents as yet unknown interactions between the symmetric TC and the specified envi- 


ronmental flow. Subject to these definitions, substituting (2.12) into (2.10) gives 
(ae CoG, #67) Cie on aoc 


SN Ne ee eee) ee) (ae) 


16 


in which the term V,+V’{, 1s absent since V, is by definition orthogonal to V’£,. Re- 
moving the terms associated with the svmmetric TC and interaction flows from (2.13) 


gives - 





( 5 = cv). ay ee) = 0, (2.14) 


which defines the evolution of the environmental winds in the absence of the TC with 


respect to the moving reference frame. Subtracting (2.14) from (2.13) gives 





a 
B44 (V5 + Ve — C)M, + (Ky — CU 
ove ie tf) = = Meet. ove, + /). (25) 


in which the (‘) symbols have been omitted for simplicitv. In all subsequent analvsis, 
independent variables will be relative to the moving reference frame. The various terms 
of (2.15) have been grouped so that the left side represents a nonlinear partial differential 
equation for the evolution of the interaction flow in response to two asymmetry- 
inducing forcing terms on the right side. The first forcing term includes: 1) distortion of 
the symmetric TC vorticity field due to large-scale horizontal windshear; and 11) steering 
of the TC by the environmental winds. The second forcing term represents the dispersive 
effect of Rossby wave radiation associated with the advection of large-scale absolute 
vorticity by the tangential wind of the symmetric TC. The TC translation velocity C also 
appears on the left side of (2.15), and represents an additional unknown. 

The effect of requiring (2.14) to hold with the TC present is that the interaction flow 
niust represent all changes to both the symmetric TC and the “basic state” environment 
due to the interaction process. However, the earlier assumption that V, may be regarded 
as steady in the moving reference frame simplifies the interpretation of (2.15) to a purely 
asvmmetric interaction flow that can cause TC propagation relative to the large-scale 
environment described by (2.14). The assumption of a steady TC is regarded as a good 
first-order approximation to lessen the considerable complexity of the TC motion prob- 
lem evident even in NDBT dynamics. However, such an assumption excludes the po- 
tential impact of TC strength and intensity changes from the scope of this research. 

As noted above. the first term on the right side of (2.15) represents both steering and 
shearing of the TC by the environmental wind. As a result. the TC translation velocity 


C incorporates both a known advection by V, as determined by (2.14), and an unknown 
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propagation velocity associated with the interaction flow V,. To separate the known 
and unknown effects, a pure shearing component of the environmental wind relative to 


the center of the symmetric TC may be defined as = 
Vir. 9,0) = V(r, 6) — V(0.0,2), (2.16) 


where V,(0,0,7) represents the instantaneous value of V, at the TC center. Note that 
V(0,0,.) provides a theoretically useful definition of environmental steering for this 


model. A propagation velocity C then mav be defined by 
Ci) = Ci) — V;(0,0,0). (2.17) 
Substituting (2.16) and (2.17) into (2.15) gives 


a a. : 
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in which the unknowns are the self-advection flow field and the associated propagation 
velocity Gi 

It should be noted that (2.14) and (2.18) represent the basis for a “propagation + 
advection” track prediction model based on self-advection that is analogous to the tra- 
ditional linear versions discussed earlier. Equation (2.14) could be numerically inte- 
grated to provide V,(0,0.7) at appropriate time intervals, while a solution to (2.18) would 
provide information on TC propagation. As discussed earlier, this research will focus on 
solving and analvzing several approximate forms of (2.18). The remainder of this chap- 
ter describes the additional assumptions and concepts on which the solutions in Chap- 


ters II] and IV are predicated. 


C. SIMPLIFYING ASSUMPTIONS 
1. Adjustment time concept 
Except for special circumstances, (2.14) indicates that the environmental 
windfield will evolve with time. The resultant temporal variability of \; and ¢, makes 
(2.18) analytically intractable for use as a theoretical propagation model without some 


sunplifving assumption regarding the time dimension in (2.15). 


18 


As discussed in Chapter I, an initially symmetric TC in a NDBT numerical 
model will propagate in response to steady large-scale forcing at a nearlv steadv velocity 
and with quasi-steady asymmetric structure after a transient adjustment period of about 
24-48 h (DM; CW; FE). The transient phase results from choosing an initial TC struc- 
ture that differs from the steady-state structure. Since TC’s are continually subjected 
to environmental forcing. it may be expected that mature TC’s in a steady or slowly 
varying large-scale windfield possess an asymmetric structure that reflects nearly com- 
plete adjustment to the asymmetric forcing of the environment. Thus, the temporal 
varibility of V. and {, in (2.18) will be eliminated by regarding the t variable in (2.18) 
as a hypothetical “adjustment time” during which a “first-guess” symmetric TC will de- 
velop a quasi-steady asymmetry appropriate to the forcing defined by (2.14) at any in- 
stant. 

In a rapidly varving large-scale environment, the TC may not be fully adjusted 
to the asymmetric forcing of the environment, such as during TC interaction with a 
rapidly moving nudlatitude trough. Without high resolution observational data to ana- 
Ivze TC propagation and horizontal asymmetries under such circumstances, a quantita- 
tive estimate for the amount of error associated with a steady-state propagation model 
is not possible. Thus, it is simply noted that the present model may be less accurate in 
such circuinstances. 

2. Matched solution approach 
a. The transition radius 
The magnitudes of the wind speed and vorticity of the symmetric TC com- 
ponent are much larger than the corresponding asymmetric and environmental compo- 
nents over a significant horizontal area in the lower and middle troposphere. Since TC 
tangential winds decay to essentially zero at large radius, a region must exist in Which 
the magnitudes of the asymmetric interaction and environmental flows are as large as 
the symmetric flow. This radius is called the “transition radius,” and will be denoted by 
R,;. The value of R,; will depend on the structure of the symmetric TC and the envi- 
ronmental windfield under consideration, and would vary with height in baroclinic situ- 
ations. Subject to some additional assumptions delineated in Chapter IV, a single 
transition radius will be assumed to exist for the barotropic model used here. 
b.  Self-advection Region 
For r < R.. it is assumed that advective terms in (2.18) not involving the 


symmetric flow will be negligible. Based on such an assumption, (2.18) simplifies to 





aA 4d Won Ul, EV eee) Vig = = Vere eV ev (Cale at 
The circular area within which (2.19) applies will be called the Self-advection Region 
since mutual advection by the asymmetric interaction flow and the symmetric TC flow 
dominate the left side of (2.19). It should be noted that the Self-advection Region as- 
sumption is equivalent to regarding the interaction flow as a linear perturbation to the 
symmetric TC flow. 
c. Dispersion Region 

For r > R,, it is assumed that self-advection terms on the left side of (2.18) 

will be negligible relative to advection terms involving the environment and either the 


symmetric TC or interaction flow. Based on this assumption, (2.18) simplifies to 


ay 
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The outer area in which (2.20) applies will be called the Dispersion Region since mutual 
advection by the interaction and symmetric TC flows, which enables barotropic vortices 
to resist dispersion, is excluded. 

d. Asymptotic assumption 

In an annular region where r = R,, the assumptions that distinguish the 
Self-advection Region and Dispersion Region are not valid. However, (2.19) and (2.20) 
are suitable for a developing a matched solution for a propagation-inducing interaction 
flow field that becomes asymptotically close to an exact solution to (2.18) as r becomes 
larger or smaller than R;. Such a matched-asvmptotic-solution approach will be used 
iInjG hapten 1a 
3. Isolating wavenumber one processes 

The diagnostic analysis by FE showed that the interaction flow field (their 
“ventilation flow’) responsible for TC propagation due to just $f has a predominantly 
azimuthal wavenumber one structure. Chan and Williams (1989) have shown that the 
interaction flow maintains an predominantly wavenumber one structure in the presence 
of horizontal wind shear in the environment, even though forcing at higher 
wavenumbers is present. The absence of significant wavenumber two and higher com- 
ponents in the interaction flow will be addressed in the stability analysis of Chapter III. 

The observation that onl a wavenumber one flow structure should cause TC 


propagation is intuitively plausible since only a wavenumber one gyre will produce 
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nonzero flow across the center of the symmetric TC (Fig. 2.2). Formally establishing 
this concept and simplifving governing equations for the Self-advection and Dispersion 
Regions is facilitated by writing (2.19) and (2.20) in component form. In a polar coor- 


dinate system centered on the symmetric TC, 


yee G.. (2.214) 
V, = u,(r,6,)r + v(r,8oe, (2.21b) 
Ci) = C() [cos(9 — a) r — sin(@—a) oe], (2.21c) 


where r and © are unit vectors in the radial and azimuthal directions respectively, and C 
and 2 are the propagation speed and direction respectively. The mathematical conven- 


tion of measuring angles counter-clockwise from east has been used. In addition, let 
a = ue oi +e 8 ,(r. 8) 1, (2.22) 


in Which conventional cartesian symbology has been used for the environmental wind 
components and unit vectors. Note that in view of the adjustment time concept, (2.22) 
contains no Gependence on time. 

a. Wavenumber one Self-advection Region equation 


Substituting (2.21a.b.c) and (2.22) into (2.19) gives 











Cs Voleean r Coe , CSs5 
= ee — fu, — Ccos\P@—a a 
ar rca A C 











CLs p Crs . Clg ‘ Cle ; 
=p == COS 6 te in sf =o as p ) cos + ts, sind. (2:25) 


(a) (b) (c) (d) 


Since both ¢],/¢r and C are independent of azimuth, the propagation term 
on the left side of (2.23) is a wavenumber one process. Only projections onto 
wavenumber one by the forcing terms on the right side of (2.23) can cause TC propa- 
gation, because projections onto other wavenumbers are orthogonal to the propagation 
term. The validity of this result depends on the assumption that the term c. Yo, may 
be neglected in both the Self-advection and Dispersion Regions. Since £, is by definition 


a function of &. projections of forcing onto other wavenumbers can potentially influence 














Vig. 2.2 Hypothetical streamfunction patterns (solid, positive; dashed, negative) for 
flows that are purely azimuthal wavenumber: (a) one; (b) two; and (c) three. Arrows in- 
dicate general direction of the flow at that location. 
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TC propagation via this term. However, FE’s finding that TC propagation was almost 
completely accounted for by advection by the wavenumber one gvre provides numerical 
support for regarding the interaction flow as a linear perturbation to the symmetric TC. 

Simphifving (2.23) to include only the dominant wavenumber one influence 
of each term 1s facilitated by writing the environmental flow factors as truncated Taylor 


series of the form 
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Since (2.16) requires A, in (2.24) to be zero for u; and t,, the forcing terms (a) and (b) 
of (2.23) project onto wavenumbers zero and two to first order. In contrast, A, is gen- 
erally nonzero for the first derivatives of environmental] vorticity. Thus, to first order the 
Self-advection Region equation governing TC propagation is 
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The parameter £,,.1s the magnitude of the environmental absolute vorticity gradient (7.¢., 


an “effective” f) defined by 
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b.  Wavenumbcr one Dispersion Region equation 

A wavenumber one propagation term does not appear explicitly in (2.20). 
Nevertheless, it may be anticipated that only wavenumber one processes in the 
Dispersion Region are important to TC propagation since orthogonality between 
wavenumbers would preclude the matching of higher order wavenumber solutions in the 
Dispersion Region to the solutions of the wavenumber one Self-advection Region 
equation (2.26). In addition, FE and Chan and Williams (1989) have shown that the 
self-advection flow is predominantly wavenumber one throughout a large region where 
the Dispersion Region assumption is applicable. 

The analysis in Section C.3.a above has already shown that the term 
V,eVf. is to first order a wavenumber two process, and by implication that the term 
V,°V(S_ + f) is also to first order a wavenumber two process. Thus, substituting 


(2:21a,b,c) and (2:22) into ( 2.20)ieives 
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To further simplify (2.29), let u, and +, be approximated by the truncated Taylor series 
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in which sin@ and cos@ terms have been used in lieu of the phase angle @ in (2.24). In 
substituting (2.30a.b) into (2.29), only terms involving sin?@ and cos’@ can contribute to 


wavenumber one. Retaining only the wavenumber one components of such terms gives 
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Noting that D, and Z, represent first-order estimates of environmental divergence and 
vorticity respectively in the vicinity of the TC, and assuming that D,x0 simplifies (2.31) 
to 





Se os . 
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which 1s to first order the Dispersion Region equation governing TC propagation. Since 
Z, 2 represents a constant angular velocity for the environmental wind, the second term 
on the left side of (2.33) represents a gvre-rotating influence consistent with the numer- 
ical results of Chan and Williams (1989). 


D. SUMMARY AND NONDIMENSIONALIZATION 

In summary, two equations have been developed that describe to first order the 
wavenumber one interaction flow associated with barotropic TC propagation. A trans- 
formation to a polar coordinate system: moving with the TC and a partitioning of the 
total flow into symmetric TC, large-scale environment and a TC-environment interaction 
flow have been emploved. A transition radius R; has been defined inside of which is a 
Self-advection Region and outside of which is a Dispersion Region where mutual 
advections by the symmetric and interaction flows are important and unimportant re- 
spectively. The wavenumber one interaction flow within the Self-advection Region 1s 


governed to first order by 


4 st Vo« Cues 
Cle vee 











ae AC. cos(6 — 0) = vsBeye Sin(O — 6), (2.34) 


and within the Dispersion Region is governed to first order by 
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where the asterisk subscripts denote dimensional variables. In all subsequent analysis, 
the absence of an asterisk will denote nondimensional variables. 
This preliminary development is concluded by deriving a nondimensional form of 


(2.34) and (2.34) by sealing r by the radius of maximum winds (R,,-), all velocities by the 
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sulting Self-advection Region equation is 
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and the Dispersion Region equation is 
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ik. BYROTROPIC VORTER STABILITY = 


In this chapter, the specific mechanism responsible for barotropic vortex stability to 
asymmetric forcing will be identified and analyzed. Since self-advection has been clearly 
identified to be necessarv for this stability to exist in barotropic numerical models (CW), 
the scope of this chapter will be limited to the Self-advection Region as governed by 
(2.32). Despite the already simplified form of (2.32), a complete solution for an arbitrary 
initial condition that includes both the transient and quasi-steadvstate responses in TC 
propagation models is in all likelihood analvucally insolvable. As a result, the focus in 
this chapter will primarily be on the transient adjustment process by transforming (2.32) 
into the related unforced initial value problem. Equation (2.32) includes only the dom- 
inant terms describing the wavenumber one interaction flow associated with barotropic 
propagation of TCs. To explain the absence of higher wavenumber asymmetries in TC 
propagation, the following development will address vortex stability to initial asymmet- 


ric perturbations of arbitrary wavenumber. 


A. MODEL DEVELOPMENT 
I. Background 

The motivation for the following analvtical development is the study of NDBT 
f-plane Couette flow by Case (1900). He obtained an integral solution for the time evo- 
lution of linear perturbations imposed as initial conditions on the NDBT Couette flow, 
which is a steady, zonally umform flow with constant latitudinal shear (Fig. 3.1a). 
Case s result mav be interpreted as an infinite summation of a continuum of singular 
solutions, hereafter referred to as continuous spectrum modes (Pedlosky 1964). For the 
NDBT f-plane Couette flow problem. this conunuous spectrum forms a complete basis, 
since discrete normal modes are eliminated by the lack of an environmental vorticity 
gradient. The superposition of these continuous spectrum modes results in a perturba- 
tion streamfunction structure that has an algebraic time dependence (i.e.. depends on 
factors involving t to an integer power) as the perturbation is ulted down-shear by the 
Couette flow. Although the baroclinic Couette flow study of Farrell (1982) showed that 
initial growth of the perturbation is possible for particular initial conditions. both Farrell 
and Case showed that the response is asvmptotically proportional to 1-? for 1 00. as 


schematically porsraved in Fag 3.Ja. In terms of instability theory, NDBT /-plane 
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Fig. 3.1 Schematic portraval of perturbation damping with time due to (a) a 
meridionally sheared Couette flow, and (b) a radially sheared axisymmetric vortical flow. 
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Couette flow may be viewed as a barotropically stable “basic state” with respect to linear 
perturbations. 

The NDBT f-plane Couette flow model, as in the Eady or Charney models, may 
be viewed as an idealization of the response of synoptic disturbances to a particular 
planetary-scale flow. Along with other limitations, the accuracy of this approach de- 
pends on the extent to which synoptic-scale disturbances may be regarded as linear 
perturbations. In the case of an intense vortex such as a tropical cyclone, such an as- 
sumpuon would be clearly unjustified since the TC winds can be substantially stronger 
than the environmental flow at quite large distances from the center. However, the large 
magnitude and nearly circular structure of the TC windfield suggests that a model anal- 
ogous to NDBT f-plane Couette flow could be developed in a polar coordinate system 
moving with the center of the TC. In such a model, the axisymmetric component of the 
TC outside the radius of maximum winds (Ry-) may be regarded as a radially sheared 
NDBT “basic state” that is time-invariant in the moving reference frame, and the asym- 
metric component of the vortex mav be regarded as a perturbation to the symmetric 
basic state. If the axisymmetric basic state has constant vorticity, theh an initial per- 
turbation can be expected to damp completely as it is tilted in the direction of the sym- 
metric radial shear (Fig. 3.16). A vortex model based on a constant vorticity basic state 
and initial perturbations is clearly a special case, and the implications of such approxi- 
mations will be addressed explicitly below. 

2. Model formulation 

The analvsis of Chapter I Section C.3 showed that shearing of symmetric TC 

vorticity by V; is to first order a wavenumber two process. It is desirable to include this 


physicallv-based process in the present analysis by rewriting (2.32) as 
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where the first two terms on the right side have been obtained from (2.23). The gov- 


erning equation for the unforced initial value problem related to (3.1) is 
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Equation (3.2) describes the evolution of an asymmetric perturbation initially imposed 
on a symmetric “basic state” vortex in a quiescent environment and on an f-plane. With 
the addition of two homogeneous boundary conditions at specified radii, (3.2) becomes 
the homogenous counterpart to (3.1) over the enclosed domain. Since these equations 
are linear with respect to the asymmetric perturbation flow, linear differential equation 
theory can be employed to show that certain relationships exist between the free- 
perturbation response and the forced-perturbation response. This issue will be specif- 
ically addressed in Section D.1 within the context of the free-perturbation behavior 
illustrated in Section B below. 

The advection of perturbation vorucity by the radially variable symmetric flow 
in term (b) of (3.2) is analogous to the shearing process that causes perturbation damp- 
ing in the Couette flow model. In addition to manifesting vortex motion through the 
advection of symmetric vorticity, term (c) permits the propagation of neutral and possi- 
bly exponentially growing discrete normal modes on the radial gradient of f,. This term 
represents a serious obstacle to an analytical solution approach since , has a strong 
radial dependence that varies with the particular vortex structure. To facilitate analysis 
of the damping process associated with term (b). term (c) will be removed by: i) re- 
quiring that +; has a Rankine radial dependence over a horizontal area bounded on the 
inside by the radius of maximum winds (Ry-); and 1i) limiting the domain of the imposed 
perturbation to that region. To facilitate specification of initial conditions in Subsection 


3 below, the modified Rankine profile (sce Anthes 1982, p. 22) will be used 
'o = ae oe (233 


The Rankine profile is the limiting value of (3.3) as Y > 1. Substituting (3.3) into (3.2) 








gives 
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where @s represents the symmetric angular wind. 

Choosing r, to be Rankine not only eliminates anv discrete normal mode com- 
ponents from solutions to (3.2), but also may have a significant impact on the eontin- 
uous spectrum response. Both these issues must be addressed before the results of this 
model can be used to interpret the numerical studies cited earlier that used non-Rankine 
vortices. The impact on the continuous spectrum can be assessed by drawing an analogy 
between (3.2) and NDBT Couette flow on a rotating sphere in the sense that the radial 
gradient of ¢, has an influence analogous to the latitudinal variation of the Coriolis pa- 
rameter. For a #-plane approximation, Kao (1955) and Boyd (1983) have shown that 8 
has no influence on the rate of continuous spectrum damping, but rather causes the 
continuous spectrum wave to retrogress proportional to f in a manner related to the 
familiar discrete mode propagation. Because the retrogression is independent of latitude 
for a constant f. the dainping rate depends only on the magnitude of the Couette flow 
shear. Since the radial gradient of ¢, is inward in the inner part of a typical vortical flow 
and decreases rapidly with increasing radius, it is anticipated that a continuous spectrum 
wave will retrogress (.e.. clockwise for a cvclonic vortex) faster at smaller radi than it 
would at larger radu. The result would be slower damping since the retrogression would 
tend to counter the tilting induced by r,. This effect should be negligible for anv modified 
Rankine vortical flow with X=1. In Section C.2, the numerical results of McCalpin 
(1987) will presented to show that significant perturbation retrogression can occur for a 
vortex that has a large and highly variable symmetric vorticity gradient. In addition, 
comparison of the results of this model with those of FE and DM will suggest that the 
symmetric vorticity gradient of a tangential wind profile that approximates a TC causes 
only minor slowing of perturbation damping due to variable retrogression. Finally. 
analvsis of McCalpin’s results will also suggest that virtually no energy from an initial 
perturbation projects onto anv discrete modes. Thus, the use of a Rankine vortex will 
be shown to be reasonable for analysis of unforced perturbation responses. 


The model may now be solved by defining a perturbation streamfunction 
4 a ¥ 7 - 
M.= hx, RS Sa (3.5a,b) 
and assigning homogeneous boundary conditions 


Ww la 0.1) = Wy 4. 6,1) — Q, (3.6a.b) 


where a= 1 and b must be chosen to confine the perturbation to the Self-advection Re- 
gion. Although the analysis of Chapter II defined the outer boundary of the Self- 
advection Region to be a specific transition radius, here a convenient choice of-b= 10 
will be used. The justification for this approach is that the Self-advection Region de- 
pends on perturbation magnitude, and the Self-advection Region in an unforced model 
can be made arbitrarily large by choosing sufficiently small perturbations. The 


azimuthal dependence of w, will also be expressed in terms of the complex Fourier series 


Wr, Ot) = Re > vr te i (3.7) 
k=1 


Substituting (3.5b) and the A” term of (3.7) into (3.4a) and integrating with respect to 


uume gives 
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Substituting (3.7) into (3.6a,b) gives 
Wat) = V(b) = 0, (3.8c,d) 


where ¢ is the radial structure of the azimuthal wavenumber & component of the initial 
asymmetric vorticity. 
Equations (3.8a-d) constitute a fully specified boundary value problem. The 


solution may be formally written in terms of a Green’s function 
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Using standard techniques (Case 1960), the Green’s function is calculated to be 
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Finally, substituting (3.9) into (3.7) gives 
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which represents a formal solution for the time evolution of the perturbation 
streamfunction. Although exact solutions to (3.11) do exist for certain initial conditions, 
the integral generally must be evaluated numerically. A simple trapezoidal scheme will 
be used. 
3. Initial condition specification 

To evaluate (3.11). an appropriate functional form and scale for ¢* must be 
specified. It is not immediately clear what combination of radial and azimuthal depend- 
encies to associate with asvmmetries generated by the interaction of a TC with its 


environmment. The right side of (3.1) is used for this purpose by writing 
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where ¢, = Ccos a and Cy = C sin » are the zonal and meridional components of the 1C 
propagation velocity respectively. Advection by C has been included from the left side 
of (3.1) since the propagation induced bv u, is also an asvmmetrv-inducing influence (cf. 
Willoughby 1988). Since asvmmetries associated with the two environmental shearing 
terms will differ only with respect to phase. which 1s unimportant in the unforced re- 
sponse analysis. only the first term will be utilized here. Applyimg similar logic to the 


tWo propagation terms in (3.12) results in the simplified form 
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To determine the spatial structure of the forcing terms in (3.13), the radial gra- 
dient of symmetric vorticity first is written in terms of (3.3). Now a Rankine profile 
AN = 1) has zero vorticity as desired to delete term (c) from (3.2). As noted in the Section 


A.2 above. a small perturbation away from (Y= 1) may also be considered as a valid 


LP) 
Ww 


model extension since the radial shear of v; would be essentially unchanged, and the as- 
sociated damping should be little affected by the small radial gradient of £, that would 


be introduced. Thus, using (3.3) with Y=1 gives - 
ee ee (3.14) 


Second, let u; be approximated by a truncated Taylor series expansion about the present 


position of the syinmetrie cyclone center 





u(y) = Sey = Ser sin @ . (3.15a) 
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where S, represents the nondimensionalized linear shear component of the environ- 
mental wind, and any zonal dependence of u; has been ignored due to the phase argu- 
ment mentioned above. 
Considering onlv term (a) of (3.13), substituting (3.3) gives 
Cd Peg 
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sin(@ — ¢), (3.16) 


Which represents the generation of a wavenumber one asymmetry from the advection 
of environmental absolute vorticity by the svmmetric vortex. Such asymmetries repre- 
sent the large gvres that have been identified in numerical models (FE) and in observa- 
tions (Chan 1986). An initial condition with equivalent spatial structure (hereafter 


referred to as a “vorucity-induced asvmmetrv”) for use in this unforced model would be 
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where a poleward gradient of absolute vorticity has been assumed ( @ = 90°), and Z,. 


represents a scale for the perturbation vorticity that remains to be specified. 


Considering only term (b) of (3.13), substituting (3.14) and (3.15a) gives 
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in which the identity 2 sin @ cos @ = sin 2@ has been used. The right side of (3.19) re- 
presents the generation of a wavenumber two asymmetry due to linear shearing of sym- 
metric vorticity by the TC-relative environmental wind. An initial condition with 
equivalent spatial structure (hereafter referred to as a “shear-induced asymmetry”) for 


USerhere 15 
0) = — c=. (3.20) 
where it is has been assumed that S, is positive corresponding to anticyclonic shear 


equatorward of the subtropical ridge. 


Considering only term (c) of (3.13), substituting (3.14) gives 





An iniual condition with equivalent spatial structure (hereafter referred to as a 
“motion-induced asymmetry”) for use here is 


“n= => k=1, (em) 
, 


where it is assumed that c, is negative. which corresponds to the generally westward 
movement of TC's relative to environmental steering documented by composite studies 
(Chan and Gray 1982: Holland 1984). In both (3.20) and (3.22), it has been assumed 
that X is slightly less than one so that the symmetric TC flow has cyclonic vorticity. 

It is inappropriate in this simple model to assign an individual perturbation 
vorticity scale Z,. to each of the above initial conditions based on the scales of the as- 
sociated forcing terms. Instead. all but one of the following model solutions will use a 
value for Z,. such that c=0.1 to facilitate analyzing the dependence of vortex stability 
on perturbation spatial structure. The radial structures of the four initial perturbations 


defiicdsabove are shownin bigs 3.2. 
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Fig. 3.2. Radial dependence of initial perturbation vorticity (Z,.¢*) for convection- 
induced (k= 2, dot; k= 1, dash), motion-induced (chaindot), shear-induced (chaindash) 
and f}-induced (solid) asymmetries. 
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B. MODEL RESULTS 
1. Perturbation damping 
Before numerically evaluating (3.11) using the initial conditions developed 


above, it is useful to obtain exact solutions by choosing initial conditions of the form 


SOS Spare (3.23) 
r 





An asymmetry of this tvxpe may be viewed as an approximation to the wavenumber & 
component of a perturbation produced by localized asymmetric convection in a TC 
evewall cloud (hereafter refered to as a “convection-induced asymmetry”). Initial con- 
dition (3.23) for k= 1 and k= 2 1s also shown in Fig. 3.2. and results in solutions to (3.11) 
of 
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Both (3.24) and (3.25) damp algebraically with time as is characteristic of con- 
tinuous spectrum mode solutions. It should be noted that the wavenuinber two pertur- 
bation damps four Umes as rapidly as the wavenumber one perturbation. which suggests 
that the damping rate is proportional to the square of perturbation wavenumber. This 
result is analogous to the dependence of damping on the perturbation zonal 
wavenumber for plane Couette flow (Case 1960). The evolution of (3.24) and (3.25) over 


a 4 


a 2 hour period are shown in Lig. 3.3a-c and Fig. 3.3d-f respectively. Unless otherwise 
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Fig. 3.3 (a)-(c) Perturbation streamfunction (solid, positive; dashed, negative) fora k= 1 
convection-induced asymmetry at t=0, 1 and 2 hours respectively. Contour interval is 
9.6 x 10? n2?/s. (d)-(f) Same as (a)-(c) except for a k=2 convection-induced asymmetry, 
and a contour interval of 4.3 x 10? 2/s. 
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noted, all illustrated solutions use the following parameter specifications: a= 1, b= 10, 
Vy-=30m s. Ry.= 50 kmand e=0.1. The contour interval is the same for each row of 
panels to show clearly the damping process. This comparison emphasizes the-strong 
dependence of the damping process on perturbation wavenumber, and also clearly illus- 
trates how the radial variability of y, ults the perturbations downshear. 

The perturbation streamfunctions for the motion-induced and shear-induced 
asvmmetries (not shown) undergo a similar, but slower shearing and damping process. 
A quantitative comparison of damping rates is facilitated by Fig. 3.4a, which shows the 
decay of perturbation kinetic energy with time for each of the four cases. Although 
perturbations of the same wavenumber damp at nearly the same rate initially, the 
damping rates quickly diverge. As will be shown in Section C.1, this property is due 
significant differences in the radial dependences of the initial conditions, e.g., r-® for a 
wavenumber two convection-induced asymmetry as compared to r-? for the shear- 
induced asymmetry. Thus. NDBT vortex stabilitv to asymmetric perturbations is quite 
Sensitive to the radial as well as azimuthal structure of the perturbation. 

The kinetic energy decay for a vorticity-induced asymmetry is shown in Fig. 
3.45. Perturbation kinetic energy decavs by approximately $0%o in 24 h, which agrees 
quite well with the adjustinent period observed in the NDBT numerical models of TC 
motion cited earlier. The much slower damping rate in Fig. 3.4b is consistent with the 
discussion in the preceding paragraph since the perturbation vorticity associated with a 
Vorticity-induced asymmetry decreases significantly slower with increasing radius than 
anv of the pret 1ous initial conditions (Fig. 3.2). The evolution of both the total and 
perturbation streamfunction components over a 36-hour period is shown in Fig. 3.5. 
The value of ¢ has been increased to 0.5 in Fig. 3.5 to emphasize how an inttually per- 
turbed NDBT vortex is restored to axisvinmetry. By § h, the inner part of the vortex 
has already regained an essenually axisymmetric structure (Fig. 3.5b), while the outer 
vorteX remains appreciably distorted. No analogy to this behavior exists in horizontal 
plune Couette flow, since the linear shear in that model renders the damping process 
independent of latitude. 

2. Influence of boundary conditions 

Boundary conditions (3.6a.bj were chosen to facilitate development of the model 
and are clearly nonphysical. Thus. it is important to determine the extent to which the 
boundary conditions may have influenced the results presented above. Two aspects of 


this prebiens must be addressed 
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Fig. 3.4 (a) Decay of perturbation kinetic energy (normalized by initial value) with time 
for convection-induced (k= 2, dot; k= 1, dash), motion-induced (chaindot) and shear- 
induced (chaindash) asymmetries. (b) Same as (a), except for a B-induced (solid) and a 
modified f-induced (dash) asymmetry. 
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Fig. 3.5 (a)-(c) Perturbation plus symmetric streamfunction for a f-induced asymmetry 
at t=0, § and 36 hours respectively. Contour interval is 7.4 x 10° »/s. (d)-(f) Same as 
(a)-(c) except showing just the perturbation streamfunction (solid, positive; dashed, 
negative), and using a contour interval of 1.7.x 108 n?/s. 
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First, boundary influences will become increasingly significant with time if there 
is a tendency for perturbation energy to propagate radially and cause a concentration 


at either boundary. This problem may be formally addressed by considering - 
C4(r, 8,1) = Reftk(rje™@-&s , (3.26) 


which is obtained by substituting (3.5b) and (3.8a,b) into the Laplacian of (3.7). 
Equation (3.26) indicates that perturbation vorticity is conserved following the symmet- 
ric angular wind, which is a purely azimuthal motion. This result depends on the model 
being linear and the symmetric vortex being Rankine, but does not depend on boundary 
conditions. Thus, no mechanism exists in this model to propagate perturbation energy 
radially, nor does an examination of Fig. 3.3 or 3.5 indicate that such a process 1s taking 
place. 

Second, the presence of the inner domain boundary at the location of maximum 
perturbation vorticity may produce nonphysical results. A nonzero inner boundarv was 
used because a Rankine vortex is singular at the origin. The inner boundary may be 
moved to the origin if the singularity is removed by modifving the wind profile to be 
solid body rotation at all radii less than Ry.. Recall that the choice of radial structure 
for the various initial perturbations depends on the symmetric vortex. Thus, the radial 
structure for a vorticitv-induced asvmmetry associated with a Rankine vortex altered as 


just suggested would be 


ae —r Vay | 
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Fig. 3.4b shows that the kinetic energy for such an initial perturbation decays at virtually 
the same rate as the original vorticity-induced asymmetry. 

Modifving the symmetric wind as described above causes the radial gradient of 
¢; to be singular at r=1. This singularity is equivalent to a large radial gradient of 
vorticity, which would tend to cause clockwise retrogression and somewhat slower 
damping of the continuous spectrum solution as discussed earlier. Although the present 
model cannot give a quantitatively precise solution for the response of a perturbation 
to a Rankine vortex with solid body inner rotation, the response of the model to (3.27) 
is useful for demonstrating that solutions are not unduly sensitive to the location of the 


inner boundary. 





C. MODEL INTERPRETATION 

1. The stabilization mechanism 
Within the context of fluid dynamical instability theory, the present stability of 
a linear perturbation with respect to some “basic state” may be assessed bv computing 
the domain-averaged local time tendency of perturbation energy. For the particular 
boundary conditions used in this model, the familiar NDBT result is that domain- 
averaged perturbation kinetic energy must be presently increasing (decreasing) with time 
if the perturbation tilts against (with) the horizontal shear of the basic state. Such a 
result is readily obtainable for this model by substituting (3.4b) into (3.4a), multiplying 
by y,, integrating over the domain, and performing a number of integrations by parts. 


Using (3.5a), the result may be expressed as 
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Equation (3.2Sa) is analogous to the familiar result for stationary cartesian coordinates 
(e.¢., Farrell 1987), except that in this case it is the shear of the axisvmmetric angular 
wind that controls the energy transfer process. In (3.28a), a positive correlation develops 
between the perturbation momentum flux and the symmetric angular wind shear when 
initially “upright” perturbations are tilted downshear. Since the radial shear of cw, is 
negative for a cyclone, the average perturbation momentum flux must also be negative. 
This flux represents an inward transport of cyclonic momentum that tends to accelerate 
the symmetric “basic state” at the expense of asymmetric perturbation kinetic energy. 
Because the analvsis thus far has been strictly linear, the symmetric vortex has 
been regarded as steady in the moving reference frame. However, using the azimuthal 
momentum equation and (3.28), the influence of perturbation momentum flux on the 


symmietric basic state 1s described by 
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This expression may be evaluated to lowest order in € using (3.11) and (3.26), and then 


integrated with respect to time to give 


Avs = 





G Cop 
a hr = = ro) Grp) Selo) [cos ki(ws(p) — ws(r)) — 1]dp, (3.30a) 


= Os(r) — ws(p) ) 
Avs(r,t) = vs(7,0) — vs(7,0). (3.305) 


The time evolution of Av, is illustrated in Fig. 3.6 using the same initial condition as in 
Fig. 3.5, except that ¢ has been reduced to 0.1, which corresponds to a more realistic 
asymmetric windspeed of approximately 2.8 m’s near the domain center. Although a 
nearly complete transfer of perturbation kinetic energy occurs by 48 h, the majority of 
the transfer has taken place by 24 h as anticipated from the damping rate in Fig. 3.4b. 
The domain-averaged sign of Ar, is positive as required, but a slight decrease is evident 
for r26. Interestingly, the radius at which Av, changes sign corresponds precisely with 
the maximum in perturbation streamfunction at any time (see Fig. 3.5f). It is also evi- 
dent that the energy transfer begins initially at small radu and spreads outward with 
tume. This aspect will be explained later in this subsection. 

The momentum flux asssociated with the convection-induced, motion-induced 
and shear-induced asymmetries (not shown) had a similar impact on the symmetric basic 
state, with the changes in Ar, becoming smaller, more concentrated at small radi and 
occurring over a shorter period of time for an initial condition that decreases more rap- 
idly with radius. The extent to which perturbation flux tends to alter the basic state is 
one measure of model linearity. Thus, the results here show that the assumptions made 
in Section A.3 (e < 0.1. r< 10) were reasonable. This demonstration of how perturba- 
tion flux transfers energy to the svmmetric vortex will be particularly useful in Section 
D.2 for understanding how a barotropic vortex achieves a quasi-steady asymmetric 
structure in the presence of steady asvmmietric forcing. 

The availability of exact solutions to this model permits valuable insights into 
the local dynamics of the stability process that are not readily evident from the 


“domain-averaged” analvsis conducted above. Note that by virtue of the identity 
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lig. 3.6 Change in svmmetric vortex windspeed with time (see legend) as a function of 
radius due to a convergence of momentum flux associated with a B-induced asymmetry. 
The windspeed of the mean symmetric vortex is 50 m/s at the radius of maximum winds 
and decreases with inverse dependence on radius to 5 m/s at ten times the radius of 


maximum winds. 
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the first and third terms on the night side of (3.24) and (3.25) have no vorticity, and thus 
are necessary solely to satisfy the boundary conditions. Therefore, substituting either of 


the exact solutions into (3.5b) results in a form 


ite V Reape s'), (3.32) 


where A(r,f) is an appropriate amplitude function based on the second term on the right 
side of either (3.24) or (3.25). Since w, is a function of radius, (3.32) will include an ex- 
plicit time dependence in some of the terms associated with radial derivatives of the 
Laplacian. For sufficiently large t, only the term possessing the highest order explicit 


dependence on t needs to be retained, which gives 
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in which the term generated by the azimuthal derivatives in the Laplacian has been re- 
tained for comparison. Note that the expression k r inside the brackets represents the 
inverse of the perturbation azimuthal length scale at radius r. If the first term in 


brackets is treated similarly, and the length scales are defined as 
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The local dynamics of the damping process may be explained in terms of (3.35) 
as follows. Under the influence of the radial shear of w;, the radial length scale of the 
perturbation at any point in the domain exhibits an inverse dependence on t, while the 
azimuthal length scale remains unaltered. Since this model conserves perturbation 
vorticity, the left side of (3.35) maintains a constant magnitude. A continuing decrease 
in L, thus requires the amphtude of the perturbation streamfunction to decrease pro- 
portional to r? so that the right side of (3.35) maintains a constant magnitude. This 
constraint may also be argued conceptually from the viewpoint that vorticity is a veloc- 


itv change over some length scale. [If the length scale is decreasing. then the velocity 
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change niust also decrease to conserve vorticity. This process is shown in Fig. 3.7 which 
is a ume sequence of perturbation vorticity using (3.26) with a vorticity-induced initial 
asymmetry. The reduction in the radial length scale due to the shearing process is clearly 
illustrated by the decreasing radial spacing of the isolines of vorticity with time. 

The results in Section B.! show that vortex stability to asymmetric perturba- 
uons is strongly dependent on the spatial structure of the initial perturbation. From 
(3.34b), this behavior can be associated with two factors that influence L,’s inverse de- 
pendence ont. The first is the perturbation wavenumber k, which results in perturba- 
uons with a higher wavenumber damping faster than those with a low wavenumber, if 
the radial structures are similar. This is the case for the wavenumber one and two 
convection-induced asymmetries (Fig. 3.3), which differ little in radial dependence (Fig. 
3.2). The second factor in (3.34b) is the radial shear of the symmetric angular wind, 
which is strongly dependent on radius. Thus, the shearing process reduces Lz, much 
more rapidly in the inner part of the vortex relative to the outer regions, which explains 
why the perturbed vortex in Tig. 3.5b has achieved an essentially axisymmetric state in 
the inner region while the outer region remains distorted. This shearing distribution also 
explains why the speed of the damping process depends on the radial structure of the 
perturbation vorticity. An initial perturbation that decreases in magnitude more slowly 
with increasing radius has a greater fraction of its kinetic energy at larger radi than does 
a perturbation that decieases rapidly with radius. Since the radial shear of w, decreases 
rapidly with radius, it will take longer to transfer the same amount of kinetic energy from 
the first perturbation. Nevertheless, energy transfer still takes place rapidly near the in- 
ner boundary. This is illustrated in Fig. 3.6 since the symmetric vortex windspeed in- 
creases alinost immediately near the inner boundary and is followed by speed increases 
at larger radii with increasing time. 

2. Verification with independent numerical results 

Insight into the impact of using a Rankine basic state vortex in this model can 
be gained by comparing the present results with a similar model that uses a non-Rankine 
vortex. McCalpin’s (1987) quasi-geostrophic reduced-gravity numerical model of an in- 
itially perturbed ocean eddy will be used for this purpose. The difference between the 
dynamical frameworks of the two models is not a significant issue since this NDBT 
model can be readily reformulated as a quasi-geostrophic reduced-gravity model. The 
result is that {, is replaced by perturbation potential vorticity in (3.4a) and that vy, be- 


comes a Kk. exponential Bessel function having zcro potential vorticity. 
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Fig. 3.7 (a)-(c) Perturbation vorticity (solid, cyclonic; dashed, anticyclonic) for a 
f-induced asymmetry at t=0, | and 2 hours respectively. Contour interval is 
pig bo KURA S 


The main difficulty in comparing the two models is McCalpin’s choice of 
Gaussian radia] dependence for both the symmetric and asymmetric components of the 
eddy. Additionally. McCalpin characterized the damping of perturbation energy in 
terms of an exponential decay time scale, whereas the continuous spectrum response in 
this model is algebraic in time. Nevertheless, if the convection-induced asymmetries used 
here are considered to be roughly equivalent to Gaussian perturbations, then 
McCalpin’s wavenuinber two decay time scale of 1.5 times the symmetric flow circu- 
lation time (at the radius of maximum Velocity) may be compared to the 0.8 h “expo- 
nential” time scale given by Fig. 3.4a for a wavenumber two convection-induced 
asymmetry. Using the 1.75 h circulation time for v,. at r= Ry. in this model, the 
damping rate here 1s approximately 3 times faster than in McCalpin’s model. This dif- 
ference is not surprising since a Gaussian vortex has an extremely strong and rapidly 
decreasing radial gradient of symmetric vorticity just outside R,,, that should cause a 
significant, radially variable retrogression that can substantially reduce the rate of per- 
turbation ulting by the symmetric flow. This assertion is supported by McCalpin’s ob- 
Servation that perturbations were advected around the eddy at only 20% of the 
maximum tangential velocity. Tangential winds in a TC decay much more slowly with 
increasing radius than a Gaussian vortex, and thus will have a symmetric vorticity gra- 
dient that is comparatively much smaller and decreases more slowly with radius. Thus, 
it mav be reasonably argued that the damping rates of this model are more represen- 
tative of those associated with typical TC wind profiles. This argument is also suported 
by the good agreement between the 24 h adjustment period cited in NDBT numerical 
studies of TC motion on a f#-plane and the kinetic energy transfer rates in Figs. 3.4b and 
3.6. With regard to the influence of perturbation wavenumber, McCalpin’s wavenumber 
three decay time scale is 2.6 times shorter than the wavenumber two time scale, which 
indicates a wavenumber dependence similar to that above. 

Earlier it was noted that choosing Rankine symmetric vortex excludes poten- 
tially important discrete normal modes. When McCalpin’s model was run on an /-plane, 
approximately 99%» of the energy initially present in the imposed perturbation was 
transferred to the symmetric component. Such a response indicates that virtually no 
perturbation energy was projected onto any discrete normal modes that might exist due 
to the presence of a symmetric vorticity gradient. Thus, the exclusion of discrete mode 
processes by the choice of a Rankine vortex for this model appears justifiable for tran- 
sient perturbation analysis. Ilowever, this does not rule out the potential unportance 


of discrete modes in the steady-state component of a forced perturbation response. 
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D. MODEL APPLICATION 
1. Free versus forced system relationships 

The motivation for the present modeling approach has been to exploit the ana- 
lytical tractability of an unforced model for analyzing transient responses. The closed- 
form solutions obtained via this technique have permitted a rigorous and illuminating 
analvsis of the perturbation damping process. To apply these results to the forced 
problem, it 1s necessary to establish what aspects of a forced-perturbation response can 
be reasonably inferred from the free-perturbation response. 

First, 11 can be shown by superposition that the full solution of a linear partial 
differential equation with steady forcing and nonhomogeneous boundary conditions can 
be represented as a sum of a steady-state part that satisfies the original equation with 
forcing and boundarv conditions, and a transient part that satisfies the homogeous 
equation, and that results from an initial condition that is different from the steady-state 
solution. Implicit in such a partitioning of the full solution is the assumption that the 
dynamical system actually supports a non-trivial steady-state condition. In the case of 
Vorticitv-induced asymmetries, the NDBT numerical study of FE indicates that the 
asymmetric component of the vortex does tend toward a steady-state condition over the 
region where the symmetric vortex is significantly stronger than the asymmetric com- 
ponent (1.e., over the region where the present linear model is valid). Although DM did 
not comment on vortex asymmetries, his illustration of a slowly varymg TC motion 
track after 24 h indicates that vortex asymmetries associated with steady, spatially vari- 
able environmental winds also tend toward a quasi-steady-state. Thus, the properties 
(e.g., damping time scales, azimuthal wavenumber dependence, etc.) of the unforced 
transient responses shown here may be expected to be relevant to the initial adjustment 
phase of a symmetric NDBT vortex on a f-plane with steady environmental winds. 

Second. in a linear dynamical system that is first order in time and exhibits a 
damped transient response to steady forcing, the magnitude of the steady-state condition 
can be expected to be proportional to the magnitude of the forcing, but inversely pro- 
portional to the magnitude of any parameter that acts to increase the rate of transient 
damping. Such an assertion 1s formally verifiable in the case of a constant coefficient 
ordinary differential equation, and may be reasonably extended to NDBT vortex dy- 
nannics for situations where numerical evidence (i.e.. FE; DM) confirms a damped tran- 
sient response to steady asvmmetric forcing. Thus, the impact of perturbation azimuthal 
and radial structure on the damping process in this model can be expected to influence 


the steady-state response to asymmetric forcing. 
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2. Barotropic vortex adjustment to steady forcing 

The process by which vortex adjustment to steady asymmetric forcing occurs in 
the various NDBT simulations of TC motion (Anthes and Hoke 1975; Kitade 1981; DM; 
CW; FE) may be explained as follows. As the initially symmetric vortex advects plane- 
tary and environmental vorticity or is distorted by environmental wind shear and vortex 
motion, an asymmetric component of vortex vorticity is generated from the symmetric 
component as indicated in (3.12). This process is a transfer of kinetic energy from the 
synimetric vortex to the growing asymmetry. If the self-advection process is omitted as 
in CW, then the energy transfer continues unchecked and causes rapid dispersion of the 
vortex. However. the analysis in Section C.1 showed that advection of the asymmetric 
vortex vorticity by the radially sheared symmetric angular wind acts to transfer pertur- 
bation energy to the symmetric vortex (Fig. 3.6) in this model. In the NDBT models just 
cited, the rate of energy transfer to the symmetric vortex apparently grows as the forced 
asymmetry grows until a quasi-steady balance is achieved over the region in which the 
symmetric vorticity is significantly greater than the asynimetric vorticity. This balance 
is in large measure achieved by about 24 h for asymmetries associated with vortex 
advection of environmental absolute vorticity (Fig. 3.4b), and occurs much faster for 
asymmetries associated with distortion of the vortex by the environment, vortex motion, 
Or as\mmietric convection (Fig. 3.4a). This variability in adjustment time scale for dif- 
ferent asymmetries is due primarily to the difference in radial] dependence between 1, 
(7!) and the radial gradient of ¢, (2-3) from which the asymmetries are generated. Al- 
though these particular dependences are specific to a near-Rankine vortex, the principal 
apphes generally since the symmetric vorticity gradient must decrease faster than the 
svmimetric Wind for all vortical flows that tend toward zero with increasing radius. 

The steady-state phase and amplitude toward which the asymmetric structure 
tends clearly cannot be addressed with a homogeneous model that also excludes the 
potentially important influence of discrete normal modes. However, two aspects of the 
steady state may reasonably be inferred. First, the steady-state asymmetry will likely re- 
tain a down-shear tilt to maintain the vortex against the continuous dispersive effect of 
the asymmetric forcing. Such a feature was noted by FE in the structure of a vorticity- 
induced astimmetry. The second inference concerns the combination of radial and 
azimuthal dependence that is likely to be present in the steady-state asymmetry. In the 
absence of environmental winds, the fi-effect results in a quasi-steady vortex asymmetry 
that is essentially wavenumber | in structure (IL). However, as shown in Chapter I], 


a horizontally variable environmental windfield will act to induce higher wavenumbers 


through distortion of the symmetric vortex. The spatial structure of the resultant 
asymmetry will depend on both the spatial structure of the environmental forcing and 
the dependence of the damping mechanism on perturbation structure. The author is not 
aware of anv detailed analyses of the wavenumber distribution and radial structure of 
TC asymmetries in either composite observations or model runs using realistic 
windfields. Thus, it is merely noted that the wavenumber and radial dependences of the 
stabilization mechanism shown here should contribute to the predominance of low 


wavenumbers and increasing axisymmetry toward the vortex center respectively. 


E. SUMMARY AND DISCUSSION 

The principal objective of this chapter has been to identify the asymmetry-damping 
influence of symmetric angular windshear as the mechanism by which a NDBT vortex 
counters dispersive and distorting influences over the region dominated by nonlinear 
self-advection. The present “linear” model captures the essence of the self-advection 
process by linearizing with respect to a nonzero symmetric basic state. In the NDBT 
models cited in Section D.2, the asymmetry-damping mechanism acts as a negative 
feedback process in which a kinetic energy transfer from the asymmetric to the sym- 
metric component of the vortex occurs as a result of, but in opposition to, the kinetic 
energy transfer from symmetric to asymmetric component induced by external forcing. 

In the previous section, the stabilization mechanism was applied to the initial ad- 
justment of a symmetric NDBT model vortex subjected to steady asymmetric forcing. 
However, the principle can also explain NDBT vortex adjustment to changes in the ex- 
ternal forcing with time. For example, assume that a quasi-steadv vortex asvnimetry 
exists due to previously steady asvmmetric forcing, and that a change now occurs in the 
environmental windfield. e.g., in f,,, (3.16) or S; (3.19). If the change is such that the 
magnitude of asymmetric forcing at a particular wavenumber is reduced (increased), then 
the shear-induced feedback of energy from the existing asymmetry to the symmetric 
vortex will be greater (less) than the environmentally-forced transfer of energy from the 
symmetric vortex to the asymmetry. As a result, the vortex will adjust toward a less 
(more) asymmetric state at that particular wavenumber until a quasi-steady balance is 
reestablished. A similar adjustment process would take place if the radial distribution 
of the forcing at anv wavenumber is altered by changes in symmetric vortex structure. 
If the duration of the external forcing change 1s brief compared to the time scale of the 
stabilization mechanism (/.e., approximating a step-function), then the adjustment time 


should be on the order of the stabilization time scale. Converselv, if the forcing is slowly 
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Varying in time compared to the stabilization time scale (e.g., a TC propagating 
poleward through a steady, but latitudinally variable windfield as in DM), then the ad- 
Justment process should have a time scale appropriate to the “apparent variability” of 
the environment from a reference frame moving with the vortex. Such a scenario is not 
intended to be all inclusive, since dynamical situations may exit in which the vortex 
might be barotropically unstable to asymmetric forcing, or in which temporary contin- 
uous spectrum growth might occur analogous to the temporary baroclinic growth 
mechanism studied by Farrell (1982). 

Symmetric angular windshear outside the radius of maximum winds can be expected 
to exert a dominant influence on the stability of any barotropic vortex. Shapiro and 
Oovama (1989) have shown that divergence had negligible influence on TC propagation, 
and thus stability to asymmetric forcing within a barotropic context. In baroclinic 
model vortices or in a TC, the role of a barotropic stability mechanism in influencing 
Vortex asymmetries will depend on the competing influence of the inertial 
stability instability made possible by the introduction of a secondary circulation into the 
dynamics. \Nfodeling studies of ocean eddies indicate that coupling between vertical 
modes can also be expected to alter vortex stability and associated motion (McWilliams 
and Flierl 1979). Significant radial shear in the tangential winds exists to large heights 
in a mature TC (e.g.. Hawkins and Imbembo 1976, their Fig. 13: Frank 1977, his Fig. 
9). Thus, the essential element that enables the barotropic vortex stability mechanism 
to operate is certainly present in TC’s. As mentioned earlier, the qualitative similarity 
of 1C motion tracks in baroclinic models (Madala and Piacsek 1975; Kitade 1980) to 
barotropic results provides at least circumstantial evidence to suggest that baroclinic 
Vortex stability is a modification to. rather than being fundamentally different from, 
barotropic vortex stability. Since, the magnitude of TC tangential windshear decreases 
with height above the boundary layer and with increasing radius outside the radius of 
maximum winds, a barotropic stability mechanism should be most influential in the re- 


gion where the TC’s convective forcing is initiated. 
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IV. BAROTROPIC VORTEX SELF-ADVECTION 


In this chapter, an analytical NDBT model of TC propagation due to planetary and 
environmental influences will be developed from the Self-advection Region and 
Dispersion Region equations of Chapter I]. The three principal obstacles that must be 
overcome are: i) the temporal dependence of the equations; ii) the strong radial vari- 
ability of the symmetric flow variables; and iii) the additional unknowns represented by 
the speed and direction of TC propagation. 

The problem of temporal variability will be eliminated by secking only a steady-state 
solution for TC propagation. Thus, it will be assumed that the model TC is fully ad- 
justed to the asymmetric forcing of the environment at any time. The complex radial 
variability of the symmetric flow will be simplified by approximating the tangential wind 
by “piecewise-defined” function. The presence of the propagation velocity in the Self- 
advection Region equation will be addressed in two ways: i) using externally generated 
values for C and o; and it) devising an internal closure scheme that will enable the pres- 
ent model to predict C and «. 


A. MODEL DEVELOPMENT 
!. Solution for the Self-advection Region 
The analysis of Chapter I1] has shown that the damping influence of symmetric 
angular windshear will permit (2.36) to evolve toward a steady-state from an unbalanced 
initial condition. Thus, the analysis here will focus on solving 
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The propagation term appears as a forcing process in (4.1), and thus introduces the un- 
known parameters C and « as noted above. The streamfunction and complex Fourier 
series definitions of Chapter III will also be used here. Thus, substituting (3.5a,b) and 
(3.7) for k= 1 into (4.1) gives 
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in Which the identities 
cos@ = Refe"} sin@ = Re{—ie"}, (4.34,6) 


have been used. and the superscript on Y’ has been omitted since all subsequent analysis 


will deal exclusively with wavenumber one asymmetries. Equation (4.2) is put in stan- 
re nv oe. 
dard form by dividing through by —, which gives 
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To obtain closed-form analytical solutions to (4.4), a convenient but realistic 


functional form must be found for the symmetric flow. In particular, requiring 
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has the desirable property of making the left side of (4.4) equidimensional similar to 


(3.8a). It mav be readily verified that 
rs") = Ayr? + Bar", (4.6) 


has an associated vorticity gradient 
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and that the product of (4.6) and (4.7) satisfies (4.5). Thus, substituting (4.6) and (4.7) 


into (4.4) gives 
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It is important to note that the radial gradient of the symmetric component of 
the vorticity for any finite vortex such as a TC must change sign at least once. Subse- 
quent calculations of R,; will show that a vorticity gradient sign change occurs within the 
Self-advection Region for TC wind profiles typically used in barotropic models. In 
contrast, the vorticity gradient given by (4.7) is monodirectional for a single value of X, 
in a region where y, is purely cyclonic. Thus, at least two segments of v, as defined by 
(4.6) for different parameters (i.e., n=1,2, . .) must be combined to approximate a TC 
symmetric windfield. The two linear (A, and B,) and one nonlinear (X,) degrees of 
freedom in (4.6) permit such a “piecewise-defined” vortex to have continuous values of 
velocity and vorticity, but a only piecewise-continous vorticity gradient. As a result, 
(4.8a.b) must be solved in at least two annular regions for particular values of X,, and 
then a total solution must be constructed by applying matching conditions at the 
interface(s) as shown in Subsections A.3 and A.4 below. The approach here is to cir- 
cumvent the difficulty in obtaining closed-form analytic solutions (¢.e., all derivatives 
continuous) to (4.4), by seeking piecewise-analytic solutions based on a piecewise- 
analytic approximation for the symmetric TC. 

Consider one of the annular regions with r=a and r=b representing the inner 
and outer boundary respectively. Within this region, a solution of (4.8a) may be ob- 
tained by noting that the functions r=*%.-' are integrating factors for the left side. 


Multiplying (4.8a) by r*.-! and integrating from a to r gives 
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and multiplying (4.8) by 7-*»~! and integrating from r to b gives 
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in which a n-subscript has been added to f to denote the nm” annulus. Performing the 


indicated integrations and eliminating ¢'Y,/¢r from the equations gives 
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in which G, results from the right sides of (4.9) and (4.10), and is defined by 
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Solutions for as many annulae as necessary to approximate the symmetric wind profile 
with segments defined by (4.6) mav be constructed from (4.11). 
2. Solution for the Dispersion Region 
Since (2.37) is a first-order hyperbolic equation in variable €,, the method of 


characteristics may be used to give an initial value solution of 
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in which {,(r. 0.0) = 0 has been assumed. The response is an undamped oscillation with 


period 4z Z, about the steady-state solution to (2.37), which 1s 
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It should be noted that no steady solution exists in the limit as Z,; > 0, and that (4.12) 


then reduces to 
Call, 0.8) = EVs Beg sin(@ = od) . (4.14) 


with £..= 6 and @= 90". 

In view of abundant numerical evidence for quasi-steady TC propagation, the 
physical relevance of (4.12) and (4.14) for large t is doubtful. Such continually evolving 
solutions are consistent with the Dispersion Region assumption that completely excludes 


the stabilizing influence of self-advective processes. Within the present modeling con- 
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text, asymmetry damping due to self-advection can impact the Dispersion Region only 
via dynamical matching conditions at the transition radius. Consequently, the evolution 
of (4.12) will limited by determining an “adjustment time” (T,,,) such that the Dispersion 
Region solution obtained below and the solution for the outermost annulus in the Self- 
advection Region satisfies an appropriate matching condition at the transition radius 
(Subsection 3). 

Making the substitution t=T,,, and applying the streamfunction (3.5a,b) and 


complex Fourier series definitions (3.7) to (4.12) gives 
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Given an appropriate form for vy, in the Dispersion Region, (4.15) may be solved using 


the integrating factors r? and r°. A convenient choice is 
26 Tp , (4.16) 


since A, and B, provide the necessary degrees of freedom to match both the value of the 
function and the first derivative of (4.6) and (4.16) atr=R,. Also, the choice of expo- 
nents in (4.16) ensures closed-form solutions to (4.15) that remain bounded as r— co. 


By the same procedure used for (4.11), the solution to (4.15) is 
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Only one annulus will be used to represent the Dispersion Region. Thus, r=a is the 
transition radius and r=b is the outer boundary of the model in (4.17a,b). In the case 


of a quiescent environment, evaluating the limit of (4.17b) as Z, > 0 gives 
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3. Matching annular solutions 
The solutions presented in Section B below will use at most three annular re- 
gions to define the model domain. Fig. 4.1 gives a schematic view of the model domain 
and annular subregions, and also shows the notation to be used to denote radial 
boundaries. Limiting values for R, and R, are the origin and infinity respectively. With 
R,=R,. two of the annulae are within the Self-advection Region, and the transition 
radius is a matching interface. The alternate use of the symbol R, for the transition ra- 
dius will permit a concise statement of the matching conditions in the following sub- 
section. 
a. Transition radius specification 
A scale analysis of (2.36) will be used to provide an empirical formula for 
R,. Implicit in the transition radius concept is the existence of an annular region about 
R; where interaction and symmetric flow variables are roughly equal. For r=R,, let 
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Substituting these relationships into (2.36) and assuming a quasi-steady balance gives 
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where the scales are shown below the respective terms. Assuming that the two self- 


advection terms do not cancel, it 1s necessary that 
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Fig. 4.1 An illustration of the model domain and three annular subdomains. The inner 
and outer boundaries of the model domain are denoted by R, and R; respectively. The 
interface between the inner (n= 1) and outer (n= 2) annulus of the Self-advection Region 
is denoted by R,. The interface between the Dispersion Region annulus (n= 3) and the 
inner annulae is denoted by R, and corresponds to the transition radius. 
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from which the empirical formula 
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(4.19) 


is obtained. Given an observed or specified radial profile of TC tangential winds and the 
magnitude of the environmental absolute vorticity gradient in the vicinity of the TC, a 
transition radius may be calculated from (4.19). 

The parameter 6 has been included in (4.19) to reflect the approximation in 
(4.18), as well as the uncertainty that will accompany any estimate of f,,. Thus, 6 is a 
“tunable” factor that could be adjusted in an “propagation + advection” forecasting 
model to minimize forecast error in some statistical sense. The extent to which a sta- 
tistically determined 6 approximates unitv would tend to provide operational confirma- 
tion of the dynamical validity of the transition radius and Self-advection Dispersion 
Region concepts. It will be shown below that assuming 6= J provides realistic results, 
and will be altered only to test model sensitivity in Section C.3 below. 

b. Boundary| Interface Condition Specifications 

Eight boundary and interface conditions must be specified to construct a 
solution for the entire model domain from solutions to (4.1 1a) or (4.1 7a) in three annular 
subregions. The asymmetric streamfunction will be sct to be zero at the boundaries of 
the mode] domain, R, and R,. and the asymmetric streamfunction and velocity will be 
made continuous at the interface radii R, and R,. The inner boundary condition is dy- 
namically required if R, is the origin. The rationale for displacing R, from the origin is 
related to functional properties of (4.6), as will be explained later. 

These conditions may be concisely expressed in terms of the complex vari- 
able ‘Y as 
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Practical implementation of conditions (4.20a-d) involves evaluating (4.1la,b) and 
(4.17a,b) at radii R,, R,. and R, to forin a 6x6 linear system of equations in the variables 
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in which double subscripts denote a matched quantity. The linear equation system will 
be solved using a standard computer algorithm. 

Because (4.8a) and (4.15a) are second order in , the associated 
boundary, interface conditions provide only enough degrees of freedom to achieve con- 
tinuous asymmetric streamfunction and velocity at R, and R,. Thus, ¢, is free to be 
piecewise-continuous at the interface radii, which is generally the case to dynamically 
balance (4.1) given the piecewise-continuous nature of the vorticity gradient associated 
with (4.6) and (4.16). As discussed above, the parameter T,,, provides an additional de- 
gree of freedom to permit matching of asymmetric vorticity at R;. In general, it will be 
possible to match only the magnitude of ¢, at R;, since varying T,,, only influences G,, 
and because the interface condition term in (4.1]7a) is irrotational. Nevertheless, the 
additional degree of smoothness given to model solution at R; should be beneficial. To 
be consistent with this quasi-matching of ¢, at R;, parameters in (4.6) and (4.16) should 
be chosen to eliminate any discontinuity in ¢f,/€r there. Not only is it possible to do 
this, but it is also desirable to facilitate closure of the model. 

4. Outline of Solution Procedure 
A NDBT solution for the propagation-inducing wavenumber one gyre 
streamifunction associated with known symmetric TC and environmental windfields has 
been obtained for an arbitrarilv large domain R, <r < R, by matching three partial sol- 
utions that are separately valid in annular subdomains. In terms of the real 


streamfunction variable w,, the solution is 
alr. 0) = Ref¥,(re"} Rp <r<R, 2 = 1,23 (4.21) 


where ‘’, and ‘’, are defined by (4.]la,b) and ‘Y; is defined by (4.17a,b). The procedure 


to compute the model results shown in the remainder of this dissertation is: 


1. Given an analvtic function that approximates a TC symmetric wind profile. deter- 
mine R, and R; (=R,) based on where the symmetric vorticity gradient changes 
sign and where the TC winds satisfy (4.19) respectively. Also, determine values for 
Ry vand Von: 


ta 


Determine parameters A,, B, and X, by constructing from (4.6) and (4.16) a 
plecewise-analvtic approximation to the wind profile used in Step 1. Because of the 
irrational function form of (4.6), this task is accomplished by successive trials using 
interactive computer graphics. The procedure is outlined in Appendix B. 


3. Assign values to the environmental parameters f,,. @ and Z;. 





4. Determine propagation parameters C and « using either external information, or a 
closure scheme internal to the dynamics of this model. 


5. Solve the &x8 linear system of equations for the value and first derivative of the 
complex streamfunction as determined by the boundaryinterface conditions 
(4.20a-d). 


6. Iteratively compute the complete model solution (4.21) by varying T,,, until the 
magnitude of asymmetric vorticity is continuous at r= R,;. 

B. MODEL RESULTS PART I: EXTERNAL CLOSURE 
As indicated in Chapter I, many observations of TC propagation are clearly sug- 
gestive of the barotropic theories, although little or no information is available con- 
cerning the TC characteristics or the environmental windfields associated with those 
vectors. Additional problems are the unknown influence of baroclinic processes on TC 
propagation and possible random systematic errors associated with the compositing 
process. Thus, only numerical predictions of TC propagation are suitable for externally 
closing this model to the degree of accuracy that accurate interpretation and analysis of 
the results can be made. The NDBT results of Fiorino and Elsberry (FE) (1989) will be 
used since that study provides both an accurate calculation of TC propagation with 
which to close this model, and detailed illustrations of the interaction flow from which 
the accuracy of closure can be evaluated. Since FE used only a quiescent environment, 
the parameter assignments f,..= 8, ¢=90° and Z,;=0 will apply throughout this section. 
The response of ihe model to environmental forcing will be addressed within the context 

of the internal closure scheme developed in Section C below. 

1. Symmetric TC Specification 
The analytic function used by FE to approximate the symmetric windfield of a 
HiGris 


pee (4.22) 


in which the parameter b varies the strength of the TC, but not the maximum intensity. 
The windspeed and vorticity gradient given by (4.22) for FE’s “basic vortex” (b= 0.96) 
are shown in Fig. 4.2a and 4.2b respectively. According to (4.18), the function #17 
(Fig. 4.2a) intersects 1,(r) at the transition radius R;= 4.75. The sign change of vorticity 
slope in Fig. 4.2b provides a good first estimate for R,. The location of the inner 
boundary R, depends of properties of the piecewise-analytic wind profile as described 


below. 
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Fig. 4.2 (a) Radial profiles of the analytic TC windfield for the Fiorino and Elsberry 
(1989) basic vortex (dashed) and the piecewise-analvtic function used (solid) to approx- 
imate it. Parameters that define the piecewise-analytic function according to (4.6) and 
(4.16) are shown in the inset. Vertical dotted lines (left to right) correspond to the radial 
boundaries, interfaces R,, R, and R, (=R,). The chain-dashed curve represents the 
function f},,.1?, (b) The analytic and piecewise-analytic vorticity gradient profiles asso- 
ciated with the wind profiles in (a). 
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A piecewise-analytic wind profile is defined by (4.6) and (4.16) for the parame- 
ters shown in the inset in Fig. 4.2a. Because the functional properties of (4.6) do not 
permit the piecewise-analvtic wind to have negative curvatire, the two wind curves de- 
part markedly near the radius of maximum winds. As a result, the values for R, and 
vs(R,) have been chosen to improve the fit of the two curves at larger radii at the expense 
of poor fit at smali radii, which is justified based on the insensitivity of TC propagation 
on TC intensity demonstrated by FE. It should be noted that the negative curvature 
region given by (4.22) near r=1 is not generally representative of TC winds, which are 
much better fit by modified Rankine profiles just outside the TC eye (Anthes 1982; p. 
24). Adjusting the piecewise-analvtic profile to accomodate such an uncharacteristic 
wind curvature is at present justified by the need to compare this theoretical model with 
the results of FE that were based on (4.22). 

The fit of the corresponding analytic and piecewise-analytic vorticity gradient 
profiles is illustrated in Fig. 4.2b. As discussed earlier, the piecewise-analytic curve has 
a discontinuity at R, due to the change in X, in (4.7) at that location. In contrast, the 
Plecewise-analvtic vorticity gradient is nearly continuous at the transition radius 
(R;=R,) in anticipation of matching the magnitude of the wavenumber one gyre 
vorticity there. The choice of R,=3.5. instead of where the analytic vorticity gradient 
changes sign (r=3), was made to facilitate the matching of gvre vorticity at R;. As in 
Fig. 4.2a, the piecewise-analvtic curve in Fig. 4.2b accurately approximates the analvtic 
profile at larger radii at the expense of a good fit at smaller radii. Because the vorticity 
gradient depends on the second derivative of the corresponding windfield, two very 
similar wind profiles can have quite dissimilar vorticity gradients. This fact has impor- 
tant implications for the dependence of TC propagation on outer wind strength, and will 
be addressed in Section C.3.b of this chapter and in Chapter V. 

2. f-induced gyre structure 

Using an initial symmetric TC wind defined by (4.22) with b= 0.96, and choos- 
ing Vy.=35 ms. Ry-=100 km and fB.=2x 10-" nr's, FE noted a f-induced propa- 
gation velocity of 2.65 ms at 330° measured clockwise from North. Using their 
parameter specifications and the piecewise-analyuc symmetric TC parameters in Fig. 
4.2a, the present theoretical model for gyre structure may be “closed” by requiring C= 
2.65 ms and x= 120° (measured counterclockwise from Last). The resulting asymmetric 
streamfunction for an outer model boundary at 10,000 km (R;= 100) is shown in Fig. 
4.3 for the central 2490x2400 km of the domain. Tor comparison, the gvre structure 


extracted from FE’s numerical mode! after 48 h is shown in Fig. 4.4. 
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Fig. 4.3 Model-predicted wavenumber one gyre streamfunction (solid, positive; dashed, 
negative) using the piecewise-defined wind profile of Fig. 4.2a, and the parameter 
specifications: }’,,=40 m/s, Ry.= 100 km, f.= 2x 10-" nr's', @ =90, R,= 1.5, R,= 100, 


C.=2.65 m’s and «=120. The contour interval is 2x 10Sn?/s, and only the inner 
2400x2400 of the domain is shown to correspond to the illustration from FE in Fig. 4.4. 
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Fig. 4.4 Numerically predicted asymmetric streamfunction at 48 h due to f-induced dis- 
tortion of a TC wind profile initially defined as in Fig. 4.2a. Contour interval is 
2x 10%nr/s, and the distance between axis tick marks 1s 40 km (Fiorino and Elsberry 
1989). 
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The analytically and numerically generated gyres have central uniform flow re- 
gions that agree quite closely with regard to flow strength and orientation, although 
some distortion due to R,#0 is evident in Fig. 4.3. Such agreement is of fundamental 
importance for the analvtical gvres to be consistent with the propagation velocity input 
from the FE numerical model. The radial positions of the streamfunction extrema are 
also accurate. However, the magnitudes of the extrema are about 20% less in Fig. 4.3. 
This is still a rather remarkable degrce of agreement considering that the two potentially 


adjustable parameters R; and T,,, have not been adjusted to improve the fit. Thus, the 


od) 
accuracy of the analvtical gyre of Fig. 4.3 provides preliminary evidence that the con- 
cepts on which this model is based are dynamically sound. Further evidence will be 
provided in Chapter V by showing that wavenumber one gvres generated from this an- 
alvtical model can accurately initialize a numerical model so that steady TC propagation 
begins almost immediately. 

One of the differences between Figs. 4.3 and 4.4 is that the analytical gyres are 
symmetric relative to a line drawn between the extrema and antisymmetric relative to the 
zero contour line, whereas the numerical solution deviates somewhat from these condi- 
tions. The first-order wavenumber one approximations in the theoretical model and the 
cyclic east-west boundary conditions used in the numerical model are hkelv responsible 
for this dissimilarity. By retaining only the wavenumber one processes, the theoretical 


analysis has excluded the interaction process 





= ee — te, (4.23) 


which results in a westward propagation of the interaction flow on the gradient of 
planetary vorticity. Such a propagation could explain the stretching weakening of the 
cyclonic gvre and compressing strengthening of the anticyclonic gyre in Fig. 4.4. This 
process is in effect generating a wavenumber two component from the wavenumber one 
gyres. The cvchc boundary condition in the numerical model artificially enhances this 
process by enabling the opposing gvres to interact in the eastern portion of the 
anticyclonic gyre. which is displaced southward relative to the analytical counterpart in 
Fig. 4.3. However, FE noted essentially negligible impact of these wavenumber two 
processes on TC propagation out to 144 hrs, which is further confirmation that a 


barotropic theory of TC propagation necds to retain only wavenumber one processes. 


6S 


3. Influence of boundary conditions 

As noted above, the choice of a nonzero inner boundary for the model has 
perturbed the uniform-flow region of the analytical gvre (Fig. 4.3). It is desirable to 
eliminate such an unphysical disturbance by moving R, to the origin. For a wind profile 
defined by (4.6), it is not possible to require the associated vorticity to be continuous for 
all r>0 and also be nonsingular at r=0. As in a geostrophic point vortex, the velocity 
in (4.6) and the higher derivatives are singular at the origin. Moving R, to the origin can 
be justified by noting that it is the ratio of the symmetric vorticity gradient to the sym- 
metric velocity that influences the interaction flow within the Self-Advection Region 
(4.4). Since the ratio of (4.6) and (4.7) contributes to a regular singular point at the or- 
igin. the interaction flow streamfunction remains defined as r—0, even though r, as 
defined by (4.6) is singular at the origin. The predicted gvre structure for R,=0 (Fig. 4.5) 
elinunates the disturbed portion of the uniform flow region, but leaves the rest of the 
gire structure essentially unchanged from that in Fig. 4.3. Since the TC associated with 
Fig. 4.5 has essentially infinite intensity, this result is in effect the ultimate extension of 
FE’s numerical demonstration that barotropic TC propagation is virtually independent 
of TC intensity. The inner boundary of this model will be placed at the origin for all 
subsequent solutions. 

The outer boundary of the theoretical model (R,;= 100) is much larger than that 
tpicallv used by numerical models. The impact of choosing R,= 28, which is about the 
distance from the center to the corners of the 4000x4000 km domain used by FE, is 
shown in Fig. 4.6. Although the peak values of the gyre streamfunction have been re- 
duced about 8°o compared to Fig. 4.5. the uniform flow region of Fig. 4.6 is virtually 
unchanged. which confirms the domain size tests of Fiorino (1987). Increasing the outer 
boundary of the theoretical model had no significant effect on the gvre structure, which 
indicates the R,= 100 is effectively infinity as far as the theoretical model response is 


concerned. 


C. INTERNAL CLOSURE FORMULATION 

In the absence of external information on TC propagation under given environ- 
mental conditions, determining the correct wavenumber gyre structure depends on a 
closure scheme to predict the associated propagation. Such a scheme will be termed 
“internal closure” since it is based solely on inherent dynamical characteristics of the 
theory. The analvsis that follows is founded on the basic hypothesis that the “correct” 


propagation velocity will be associated with a particular gyre structure that is unique and 


69 





Y-DISTANCE (Ry.) 








-9 —-6 -3 ; 3 6 9 
X-DISTANCE (Ry) 


Fig. 4.57" Asin Fig. 4.3 texcept for R= 0) 
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As in Fig. 4.5, except for R,= 28. 


readily discernable. Since the work of FE will be used to provide both motivation for 
and validation of this approach, the closure scheme will be developed within the context 
of TC propagation in a quiescent environment, and then applied in Section D below to 
illustrate the response of the model to variable environmental winds. 
1. Preliminary Analysis 

As noted in Chapter I, FE found that the advection of the wavenumber one gyre 
by the symmetric winds of the TC produces a streamfunction tendency that nearly can- 
cels the tendency associated with the advection of planetary vorticity by the symmetric 
winds. Shapiro and Oovama (1989) have recently refined this concept further by show- 
ing that tangential windshear produces nearly uniform asymmetric absolute vorticity 
within =350 km of the TC center. As shown in Chapter II, this process is fundamentally 
a kinetic energy feedback loop in which a transfer of energy from the asymmetry to the 
symmetric TC due to symmetric angular windshear occurs as a result of, and in oppo- 
sition to, the transfer of energy from symmetric TC to the asymmetry induced by a gra- 
dient of planetary vorticity (/). 

The closure scheme developed here will implement the above concepts by seek- 
ing a particular balance among the various processes that contribute to interaction flow 


vorticity tendency. Using (4.1), the balance in the Self-Advection region must be 


Av Ay Ay ¥ 





ST eee Sh C= acos6 Sa 20, (42 
= ay — vshcos@ — u, a zy cos(0 — a) = = (4.24) 


in which a quiescent environment has been assumed, and the identity on the nght em- 
phasizes the steady-state nature of the present model. In light of Shapiro and Ooyvama’s 


observation, (4.24) will be scaled by — €§,/Cr to give 





1 4 as : 7 
| r Oo + fBcos | amc + u, — Ccos(@-—a) = 0, (4.25) 
(a) (b) (c) (d) 


which is valid at all radi except where the symmetric vorticity gradient goes to zero. 
Terms (a)-(c) of (4.25) represent contributions to the vorticity tendency of the inter- 
action flow associated with symmetric angular windshear, f-induced dispersion, and 
advection of symmetric TC vorticity by the uniform portion of this interaction flow re- 
spectively (hereafter the “Shear”, “Beta” and “Advection” terms). Term (d) represents a 


vorucity tendency that is an artifact of the movement of the model reference frame with 


the propagating TC (hereafter the “Motion” term). The quantity within the brackets in 
(4.25) is in fact the asymmetric absolute vorticity gradient, which according to Shapiro 
and Ooyama must become nearly zero inside about 350 km for a TC undergoing quasi- 
steady propagation. Thus, the cancellation of the Shear and Motion terms results in the 
Motion term being balanced solely by the Advection term, as observed by FE. 

Since all the terms on the left side of (4.25) are wavenumber one processes, a 
particular direction may be associated with the maximum vorticity tendency of each 
term at any radius. The vectors in Fig. 4.7a,b show the magnitude and direction of such 
vorticity tendency maxima at r=2 and 4 respectively associated with the wavenumber 
one gyre streamfunction (Fig. 4.5) predicted by this model using a piecewise-analvtic 
approximation to FE “basic” vortex (Fig. 4.2a) and their reported propagation velocity 
of 2.65 ms at 330° (x=120). Since (4.25) has been scaled by the symmetric vorticity 
gradient, the Shear, Beta and Advection terms appear in Fig. 4.7a.b as velocity compo- 
nents that must together balance the imposed Motion vector. Since the symmetric 
vorticity gradient changes sign in the Self-advection Region, two radii have been chosen 
to permit comparison of the balance between the inner (r= 2) and outer (r= 4) annulus 
of the Self-advection Region. 

A radial profile (Fig. 4.7c) of the amplitude and phase of the interaction flow 
vorticity (f,) corresponding to the streamfunction in Fig. 4.5 highlights distinct amphi- 
tude and phase changes between the inner and outer annulus of the Self-advection Re- 
gion. Identifving specific patterns in both this profile of interaction flow vorticity and 
the orientation of vorticity tendency vectors above is the basis for the following closure 
scheme. 

2. The Closure Schenie 

Because the Shear and Beta vector tend to oppose each other in Fig. 4.7a and 
b, the Advection vector balances most of the imposed Motion vector. Such a vector 
pattern may be regarded as roughly consistent with FE’s observation that the 
streamfunction tendencies associated with the Shear and Beta terms of (4.25) nearly 
cancel in their numerical solutions. If barotropic TC propagation is indeed precisely 
steady propagation, then the Beta and Shear vector would be expected to cancel exactly. 
The inexact cancellation in Fig. 4.7 may be attributed to inaccuracies in the theoretical 
model and or the possibility that barotropic TC propagation is more a quasi-steady 
process than a truly steady process. Assuming that the second factor is primarily re- 
sponsible for the vector pattern in Fig. 4.7, a closure scheme based on precisely steady 


propagation will be feasible if: 1) a propagation velocity can always be found sucli that 


~} 
la 


SPEED (M/S) 














-1 0 

SPEED (M/S) SPEED (M/S) 
——f® SHEAR TERM —> MOTION TERM 
SS re ee ean — > ADVECTION TERM 





a 
<< 
1.0 
180 


0.8 
90 


0.6 


0.4 


AMPLITUDE 


0.0 0.2 





RADIUS 


Fig. 4.7 (a)-(b) Vector diagrams showing direction and magnitude of ¢, tendency 
maxima that have been scaled by — 6¢,/6r at (a) r=2 and (b) r=4 respectively. The 
vectors represent the terms of (4.24) (see key above), and were computed using the 
piecewise-defined symmetric TC wind profile of Fig. 4.2a and the wavenumber one gyre 
solution of Fig. 4.5. (c) Radial profiles of the amplitude (solid) and phase (dashed) Ol <, 
corresponding to the streamfunction field of Fig. 4.5. The amplitude curve has been 
scaled by the amplitude of ¢, at r=0.2. Vertical dotted lines (left to right) correspond 


to Ay dnd Be. 
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theoretical model predicts precise cancellation of the Shear and Beta vectors; and ii) such 
a propagation velocity is a reasonable approximation to the quasi-steady result obtained 
from a dynamically equivalent numerical model. As shown in Fig. 4.8a. the propagation 
velocity C= 1.9, a= 130 satisfies the first of the above conditions. Such a closure scheme 
is clearly radius dependent since the Shear and Beta terms in Fig. 4.8b are opposite in 
phase, but are somewhat unequal in magnitude. Although the theoretically predicted 
Propagation direction differs by only 10° from the propagation direction reported by FE, 
the speed difference of 0.75 m s represents a significant error. 

Both the speed error and the radius dependent nature of the above closure 
scheme may be addressed by recalling FE’s demonstration that TC propagation and 
advection by V, in the uniform flow region are nearly, but not precisely equal. Thev 
observed vector differences between the propagation velocity of the TC and the uniform 
poruon of the interaction flow (their “Venulation flow”) that are persistently westward 
at about 0.3 ms (Fig. 4.9). This difference is readily evident quite early {12 h) in the 
integration. which suggests a linear Rossby wave-type propagation associated with the 
subsynoptic scale Fourier components of the symmetric TC. Such a hypothesis is sup- 
ported by Shapiro and Oovama’s (1989) Fourier-Bessel spectral analysis. which showed 
that the squared angular momentum spectrum for a typical TC wind radial profile has 
a distinct peak at a radial length scale of 200 km. If this radial length scale is regarded 
as approximately equivalent to the inverse of a sinusoidal wavenumber, then it mav be 
shown that a ssmmetric nondivergent Rossby wave of such a length scale will propagate 
westward at 0.4 m 5. 

The requirement that the Shear vector exactly cancels the Beta vector in the 
closure scheme is inconsistent with these numerical results that only the component of 
the Beta vector over about 0.4 ms actually contributes to distortion of the TC circu- 
lauon. That is, only the relative linear propagation speeds of the large and small scale 
Fourier components contributes to the distortion of the TC circulation. Thus, a better 
closure formulation is to choose C. and @ so that at a radial location 200 km from the 
TC center the Beta and Shear vectors are 180° out of phase and differ in magnitude by 
04 ms (Fig. 4.10a). The resulting values for C. and a for this closure scheme are 2.4 
ms and 132° respectively. which reduces the speed error to 0.25 ms. Such an error is 
considered quite acceptable in view of the typical magnitudes of total TC motion 
(advection ~ propagation). The Beta and Shear vectors in Fig. 4.10b also differ by 
about 0.4 ms, which results in Advection and Motion vectors at r=4 that closelv ap- 


proximate those atr=2. In contrast to Fig. 4.7c¢, the phase of interaction flow vorticity 
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Fig. 4.9 Vector differences between the TC propagation and an averaged interaction flow 
velocity (0-300 km) at various times (see legend) during a 12U h integration of a 
barotropic numerical model (Fiorino and Elsberry 1989). 
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Vig. 4.10 As in Fig. 4.7, except using the analytical model solution based on C.=2.4 m/s 
and o = 132. 
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is constant for 1.5<r< R; in Fig. 4.10c, which indicates that the vector patterns in Fig. 
4.10 are characteristic of most of the Self-advection Region. Such consistency of the 
vorticity tendency balance is to be expected for a TC propagating as a quasi-steady en- 
uty, and supports the dynamical validity of both the theoretical model and the modified 
closure scheme. 

It should be noted that the phase of ¢, changes abruptly for r<1.5. In addition, 
the amplitude of {, increases quite rapidly and actually becomes singular at the origin. 
Both of these properties are characteristic of an intense “inner gyre” in the vorticity field 
of the interaction flow, which must exist to balance the singular nature of the 
piecewise-analvtic symmetric flow at the origin. Thus, the inner vorticity gyre may be 
purely an artifact of this theoretical model. By contrast, numerical modeling studies 
suggest that such mner vorticity gvres result from a dynamical instability which is sup- 
ported by the barotropically unstable structure of typical TC tangential wind profiles 
(e.g., Willoughby 1988; Peng and Wiliams 1989). Similar gyre patterns have been sug- 
gested by aircraft observations of the inner wind fields of a TC, although these may be 
to some extent due to mislocating the 1C center. In this model, the inner gyres have 
no significant impact on the propagation prediction capability of the internal closure 
scheme just demonstrated. 

The usefulness of the internal closure scheme depends on the “closure” vector 
pattern (/.e., Fig. 4.10a) occurring at a well-defined and unique location in the C. and 
Parameter space. By showing that the vorticity tendency vector pattern changes dis- 


” 


tinctly in response to small adjustments of C. and away from the “closure” values, 
Figs. 4.11-4.14 verify that the closure point may be quite precisely located. Changes in 
o to 137° (Fig. 4.11) and 127° (Fig. 4.12) induce a meridional component into the Shear 
vector. In contrast, changes in C.to 2.1ms (Fig. 4.13) and to 2.7 ms (Fig. 4.14) tend 
to alter the relative magnitudes of the Beta and Shear vectors, although the processes 
are coupled to a moderate extent. The phase uniformity of £, in the Self-advection Re- 
gion is rapidly altered (Fig. 4.11c and 4.12c) when the specified propagation direction 
differs by only +5° from the correct value. Although all closures in this research are 
obtained via interactive computer graphics. the well-defined nature of the closure vector 
pattern should be amenable to computer automation. A formal proof that only one 
closure pattern exists for a particular set of TC and environment parameters has not 
been attempted. However, a careful search over a wide range of imposed propagation 


velocities has never revealed more than one closure point in C. and @ parameter space. 
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Fig. 4.11 As in Fig. 4.10, except using C.= 2.4 m/s and a= 137. 
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Fig. 4.12 As in Fig. 4.10, except using C.=2.4 m,s and o= 127. 
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Vig. 4.13 As in Fig. 4.10, except using C.=21 m/s and «= 132. 
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Fig. 4.14 As in Fig. 4.10, except using C.= 2.7 m/s and «= 132. 
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PHASE 


The closure scheme quite precisely locates a “predicted” direction of propagation 
that is 12° to the left of the actual value found by FE. The cause for most of this bias 
is that the phase of w, converges to a value of 9° at large radii (Fig. 4.15), rather than 
to zero, as would occur in a quiescent environment if the Dispersion Region equation 
was solved over the entire model domain. The mathematical basis for this property of 
the solution is that in (4.17a) the terms associated with the inner boundary condition and 
to first order the forcing term G,/2r both decay proportional to r forr>R,. Asa result, 
the symmetric angular windshear that produces the barotropically stable downshear 
phase shift in and near the Self-advection Region (Fig. 4.15; r<8) still has a small but 
measureable influence in the Dispersion Region at large radu. Since the phase shift of 
yw, at large radii mav always be computed from the model solution, the closure-predicted 
propagation direction could be corrected by this angle to increase the accuracy of the 
model for applications in which a precise direction prediction would be of particular 
importance (e.g., a propagation+advection track forecasting aid). Such directional 
corrections will not be used here. As will be shown in Chapter V, the corrections are 
not required when using this model to initialize a numerical model. 

The accuracy of the internal closure scheme developed above may be put into 
perspective by computing the 24 h forecast errors that would result solely from speed 
and direction errors of 0.25 ms and 12°. which represent the difference between the 
propagation predicted here via the internal closure scheme and that observed by FE. 
The speed error translates into =25 km of forecast error. Assuming a typical TC prop- 
agation velocity of 2.5 ms, the direction error can cause =45 km of forecast error. 
These errors are about one quarter of the tvpical 24 h official forecast error (e.g., 
Thompson et al. 1981), and thus would be of secondary importance compared to errors 
induced by poor, inadequate observations. 

3. Sensitivity testing 
a. Transition radius adjustments 
This model contains two potentially adjustable features: the location of the 
transition radius as determined by 6 in (4.19), and the requirement for a relatively small 
change in [{,| at r=R,;, which may be defined in terms of the parameter y by 
| oy (r,) | 


y= lim ——. 4.26 
rR | ( ) 
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Fig. 4.15 Radial profile of streamfunction amplitude (solid) and phase (dashed) corre- 
sponding to vorticity profile shown in Fig. 4.10c. The amplitude curve has been nor- 
malized by the maximum amplitude of the streamfunction. 
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Thus, it is important to test the sensitivity of the closure scheme to moderate adjust- 
ments of these conditions, and also to determine whether any physically unreasonable 
changes occur in the predicted wavenumber one gyre structure, 

Table 4.1 is a summary of the sensitivity tests performed, including the re- 
sultant changes in T,,,. Test 1 and (2) shows that decreasing (increasing) R; results in a 
less (more) westward direction of propagation, and a slower (faster) speed of propa- 
gation. In contrast, decreasing (increasing) y in Test 3 (4) results in a less (more) 
northward direction of propagation, and a slower (faster) speed of propagation. Since 
changes in y are accomplished by adjusting T,,. in (4.17b), Table 4.1 includes the values 
of T,,. to show that changes in T,,,. are approximately proportional to the changes in 
a which is consistent with the idea that a stronger outer gvre structure will contribute 
to stronger interaction flow across the center of the TC. In each test, the phase of w, 
at large radu was found to be similar in magnitude to that noted above (e.g., Test 3, 
11.1°; Test 4, 7.5°). In Test 4, the predicted propagation speed and direction are quite 
close to FE’s 2.65 m's at 330° (a= 120°). The wavenumber one gvre structure associated 


with the C. and a of Test 4 is shown in Fig. 4.16. If the wavenumber two component 


Table 4.1 Response of the model-predicted propagation velocity (columns 5,6) to four 
combinations of the parameters 6 (column 2) and y (column 3) as defined by (4.19) and 
(4.26) respectively. 
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Fig. 4.16 As in Fig. 4.5, except using the parameter specifications of Test 4 of Table 4.1. 
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is ignored, then the numerical result of FE (Fig. 4.4) is quite closely approximated by the 
theoretical result in Fig. 4.16. Thus, the sensitivity tests demonstrate that the theoretical 
model depends of the adjustable parameters in a well-defined and physically reasonable 
manner. 

b. Piecewise-analytic symmetric TC specification 

The piecewise-analytic TC winds profile (Fig. 4.2a) used above was con- 
structed so that the associated vorticity gradient (Fig. 4.2b) is continuous and equal to 
the analytic vorticity gradient at the transition radius. In contrast, the piecewise-analytic 
vorticity gradient is discontinuous at the interface R, and only approximates the analytic 
vorticity gradient in an average sense. The motivation for approximating the analytic 
vorticity profile in this manner may be illustrated by determining the model-predicted 
propagation for the piecewise-analytic profiles in Fig. 4.17. Note that the piecewise- 
analytic wind profile approximates the analytic wind profile more closely than in Fig. 
4.2a. In contrast, the piecewise-analvtic vorticity gradient is now discontinuous and only 
poorly approximates the analstic profile at both R, and R;. The model-predicted prop- 
agation using the parameters of the figure inset is C.= 1.75 and «= 144°, which deviates 
significantly from FE’s results of 2.65 ms and 330° (a = 120). 

The cause of this important error in model-predicted propagation may be 
determined by eliminating the processes that cannot be responsible. Since the 
piecewise-analytic wind profiles in Figs. 4.2a and 4.17a are nearly identical, the change 
in predicted propagation cannot be attributed to the Beta or Shear terms of (4.25). That 
is, the asymmetric forcing associated with a planetary vorticity gradient and the sym- 
metric angular windshear-induced stability of the TC remain essentially unchanged. The 
change in the piecewise-analvtic vorticity gradient in the Dispersion Region cannot be 
responsible since self-advection processes are excluded there. The rough agreement be- 
tween the piecewise-analvtic and analytic vorticity profiles inside R, are also quite similar 
in Figs. 4.2b and 4.17b, and FE have shown conclusively that changes in the symmetric 
TC inside 300 km have a negligible effect on propagation. Thus, it must be concluded 
that the significant reduction in the average magnitude of the piecewise-analvtic vorticity 
gradient in the outer annulus of the Self-advection Region must be responsible for the 
change in predicted propagation. Since the piecewise-analvtic vorticity gradient roughly 
brackets the analytic curve in Fig. 4.17b near R,, the underestimation of the analytic 
Vorticity gradient just inside the transition radius is inferred to be the source of the error 


in predicted propagation. 
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Fig. 4.17 (a)-(b) As in Fig. 4.2, except that the parameters (inset) of the piecewise- 
analytic TC wind structure have been changed so that the associated vorticity gradient 
is no longer continuous and equal to the analytic value at the transition radius (r= 4.75). 
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10 


This analysis suggests that the sensitivity of TC propagation on wind 
strength in a 300-800 km “critical annulus” demonstrated bv FE is actually a manifesta- 
tion of the sensitivity of the self-advection process to the radial gradient of symmetric 
TC vorticity via the Advection term of (4.25). Thus, barotropic TC propagation actually 
depends on the second derivative of the symmetric windfield in the critical annulus, 
which explains why seemingly modest changes in outer wind strength as used by FE can 
cause significant changes in propagation. The piecewise-defined nature of the present 
theoretical model tends to concentrate this dependence near the transition radius, and 
also exaggerate the dependence by enabling very similar piecewise-analytic wind profiles 
to have associated vorticity profiles that differ more severely than similar analytic pro- 
files. 


D. MODEL RESULTS PART II: INTERNAL CLOSURE 
The availability of an internal closure scheme permits this model to predict TC 
propagation for various TC wind profiles and environmental influences. Efforts are 
presently underway to compute TC propagation relative to horizontally variable envi- 
ronmental winds from numerical models (personal communication, R. T. Williams 
1989). However, no published results are available against which the above internal 
closure may be compared. Thus. it will be assumed that the closure scheme with 6= 1 
and y=1 will give results of comparable accuracy when applied with spatially variable 
environmental winds. An internal check of the likely accuracy of the prediction can be 
made by comparing the closure vector plot for dynamical consistency with the associated 
gyre structure predicted by the model for various environmental windfields. 
1. Influence of TC structure change 
Chan and Williams (1987; CW) also used (4.22) to specify the initial TC 
windfield, and varied TC intensity and strength simultaincously by adjusting the pa- 
rameter J’. with b=1. Three of their parameter specifications (Table 4.2; columns 2-4) 
will be used here. In light of the closure testing above, the piecewise-analvtic profiles 
(Figs. 4.18-4.20) have been constructed so that the analytical vorticity gradients are ac- 
curately approximated near the transition radius. Note also the increase in the transition 
radius {Table 4.2; colunmin 5) in response to larger valucs of J’... The comparison of the 
theoreticallv-predicted propagation speeds and directions with those reported by CW 
(Table 2: columns 6.7) reveals slow speed biases and left direction biases quite similar to 
those in Section C above, which indicates that the model results are quite consistent with 


regard to changes in TC structure. As shown in subsection 3.a above, tlie persistent 
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Fig. 4.18 As in Fig. 4.2, except using Chan and Williams (1987) analytic TC wind profile 
parameters of 1 ,.=20 m/s, Ry.= 100 km and b= 1.0. 
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Fig. 4.19 Asin Fig. 4.17, except for Vy,.=40 m/s. 
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speed and direction biases could be reduced by moderately adjusting the parameters 6 
and y. 

It should be noted that the theoreticallv-predicted propagation directions in 
Table 4.2 are essentially independent of changes in TC structure, whereas FE noted a 
more Westward propagation direction in response to increases in TC outer wind strength. 
The track of FE’s “large-weak vortex” matched closely that of the “basic vortex” for the 
first 24 h, but became increasingly westward as the integration proceeded. Thus, it 1s 
possible that the numerical model includes some non-steady propagation process that is 
excluded from this steady-state model. An alternate explanation is that cyclic boundary 
conditions in the numerical model may have a significant impact on the propagation 
vector associated with a large and purely cyclonic vortex that radiates substantial 
amounts of Rossby wave energy (cf. Shapiro and Ooyama 1989), Althougth domain size 
tests Were conducted by FE with the basic vortex, it does not appear that such tests were 
conducted for the large-weak vortex. The combined evidence of the theoretical results 
here and the numerical results of CW suggest that changes in TC strength result prima- 


rily in changes in propagation speed, but not propagation direction. 


Table 4.2 TC propagation velocity (columns 6.7) predicted by the theoretical model and 
the numerical model (in parentheses) of Chan and Williams (1987) for three values of 
maximum svmmietric Wind (column 2). The analytic and piecewise-analvtic curves used 
by the numerical and theoretical models respectively are shown in Figs. 4.18-4.20. 
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2. Influence of uniform environmental vorticity 


The next step in complexity from f-induced TC propagation in a quiescent en- 
vironment is to include a uniform environmental vorticity field. ie., Z;#0, Vo,;=0. Ex- 
amples of this situation would be: 1) approximating the trade winds as a nondivergent 
zonal flow with linear meridional shear; and ii) approximating the subtropical ridge 
(monsoon trough) region as a large-scale anticvclone (cyclone) that is in solid body ro- 
tation (cf. Willoughby 1988). 

Table 4.3 is a summary of the closure results for four different values of envi- 
ronmental relative vorticity using the same piccewise-analytic TC as in Case 2 of Table 
4.2. Recalling that closure occurs at C.=2.65 m s, a= 132° for Case 2 of Table 4.2, a 
cyclonic (anticyclonic) environmental vorticity induces a counterclockwise (clockwise) 
adjustment in the direction of TC propagation that is proportional to the magnitude of 
Z;-- Note that the theoretical model also predicts a greater TC propagation speed as Z;. 
progresses from anticyclonic to cyclonic, which suggests that TC propagation mav be a 
more dominant component of total TC motion in regions of large-scale cyclonic vorticity 


such as the monsoon trough. Although no numerical model confirmation exists for this 


Table 4.3. Theoreticallv-predicted propagation velocity (columns 3.4) for four values of 
uniform environmental vorticity (column 2). 
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aspect of the prediction, it is consistent with the faster propagation of westward moving 
TC’s in composite data (Fig. 1.1b). The wavenumber one gvre patterns associated with 
Cases 1 and 4 of Table 4.3 are shown in Fig. 4.21a and b respectively. Chan and 
Williams (1989) have noted gvre orientation and propagation changes in an equivalent 
numerical simulation that are consistent with the results given here. 

It is important to note that this model cannot be internally closed if Z;. is much 
less or greater than the range of values in Table 4.3. The mathematical explanation for 
this problem is that the solution for ¢, in the Dispersion Region depends inversely on 
Zr. For sufficiently large values of environmental vorticity, it becomes impossible to 
iterate to a value of T,,. such that the magnitude of ¢, is continuous at the transition 
radius. From a dynamical perspective, the assumptions V;«<V, and (,<{, become un- 
acceptably inaccurate in the outer part of the Self-advection Region. As a result, the 
present form of the theoretical model can accomodate only moderate values of large- 
scale vorticity, and should not be applied to situations involving interaction of the TC 
with an intense cyclonic system such as another TC. Such situations might be addressed 
by extending the model to include some of the higher-order wavenumber one processes 
that were excluded here (Chapter I]). 

3. Influence of environmental vorticity gradients 

To determine the model response to gradients of absolute vorticity, combina- 
tions of Z;., B... and @ are devised that approximate the horizontal variability of the 
typical large-scale environments through which 1C’s move. Figs. 4.22a and b illustrate 
parameter combinations that correspond to the regions south and north of the NW 
Pacific subtropical ridge respectively. As indicated, Z,;. is assumed to be moderately 
anticyclonic in both regions. The westward component of the environmental vorticity 
gradient reflects the influence of the intense Asian heat low throughout much of the 
troposphere during the NW Pacific typhoon season. Similarly, the southward and 
northward components reflect the influence of the equatorial trough and midlatitude 
baroclinic regions respectively. Such parameter combinations produce an absolute 
vorticity gradient that is larger and is directed more poleward north of the subtropical 
ridge compared to south of the subtropical ridge (Table 4.4; columns 3,4). Although the 
analvuc TC wind profile of Table 4.3 Case 2 is used for both cases here, an additional 
piecewise-analvtic profile (Fig. 4.23) has to be constructed to account for the increased 
magnitude of #,,. north of the subtropical ridge. 

If the phase shift of the environmental vorticity gradient is taken into account, 


then the tendency vector patterns (Figs. 4.24a,b and 4.25a,b) from which the propa- 
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Fig. 4.21 Theoretically-predicted wavenumber one streamfunction fields for (a) Case 1 
and (b) Case 4 of Table 4.3. Contour interval is 2 x 10°r7/s. 
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Fig. 4.22 Idealized planetary and environmental vorticity gradients for TC positions in 
an anticyclonic vorticity region (a) south and (b) north of the subtropical ridge in the 
lower to middle troposphere of the western North Pacific TC region. The units of the 
vorticity oradients are 10-"yi-'s-)) and the wniisol 77. atelier), 


98 





VELOCITY (V4) 


Oe eet sie) 





Ro vs{R,) a 


1.50 1.00 
3.25 0.33 0.10 -3.69*10" 5.05*10° 
4.25 0.16 2.50 —4.14*10° 6.43*10° 

1.03*10° —2.20*10° 


a a ae a ee 
RADIUS (Rye) 


Fig. 4.23 As in Fig. 4.18, except that the piecewise-defined wind parameters (see inset) 
have been recalculated to account for £,,.= 3.5x 10-7! ms” in Fig. 4.21b. 
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gation velocities in Table 4.4 were obtained are consistent with those shown earlier (e.g., 
Fig. 4.10). In each case, the Shear vector phase is opposite that of the Beta vector, 
which is now oriented 90° to the left of the direction of the environmental absolute 
vorticity gradient. The characteristic patterns in the radial profiles of interaction flow 
vorticity magnitude and phase are also evident (Figs. 4.24c and 4.25c). The orientation 
and amplitude of the wavenumber one gvre streamfunction are also dynamically con- 
sistent (Fig. 4.26), which suggests that the accuracy of the model with an environmental 
vorticity gradient included should be similar to the results shown earlier. Although 
equivalent numerical results are not available, the propagation directions predicted here 
are quite similar to the observed propagation directions for westward (pre-recurvature) 


and eastward (post-recurvature) moving TC’s (Fig. I.1b). 


Table 4.4 Theoretically-predicted propagation velocities (columns 5.6) in response to 
environmental parameter combinations (columns 2-4) representing TC locations south 
(Case 1) and north (Case 2) of the subtropical ridge during the western North Pacific 
typhoon season. 
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E. SUMMARY 

An analytical NDBT model based on the principle of nonlinear self-advection has 
been developed to predict the steady-state TC propagation and associated wavenumber 
one interaction flow induced by planetary and environmental forcing. The model has 
been made analytically tractable by dividing the highly complex and generally intractable 
TC propagation process into a number of manageable pieces via a number of reasonable 
assumptions. A piecewise-analytic solution for the wavenumber one interaction flow 


and the associated propagation velocity was then constructed from analytic solutions to 
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Fig. 4.24 (a)-(b) Vorticity tendency vector diagrams, and (c) vorticity amplitude and 
phase profiles as in Fig. 4.7, except for the mode} solution for Case | of Table 4.4. 
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Tig. 4.25 As in Fig. 4.23, except for Case 2 of Table 4. 
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Fig. 4.26 Model-predicted wavenumber one streamfunction fields as in Fig. 4.3, except 
for (a) Case 1 and (b) Case 2 of Table 4.4. Contour interval is 2 x 10%»7/s. 
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the individual pieces of the problem. The two key aspects of this piecewise-analytic 
modeling approach are to: 


e divide the model domain relative to a transition radius within which is a Self- 
advection Region where mutual advections by the symmetric TC and the asym- 
metric interaction flow are important, and outside of which is a Dispersion Region 
where such advections are considered unimportant; and 


® approximate the symmetric TC windfield by a piecewise-analvtic modified Rankine 
profile that closely matches the analytic TC wind profile, except near the radius of 
maximum winds, and closely matches the analytic vorticity gradient in the Self- 
advection region near the transition radius. 
This theoretical model is distinguished from previous efforts by the capabilities to: 


e accurately predict both the zonal and meridional components of TC propagation 
by integrating linear and nonlinear mechanisms into a single self-advection process: 


e accurately predict the wavenumber one gvre structure responsible for TC propa- 
gation based on either a knowledge of the propagation velocity, or determining the 
propagation from the model dynamics via a closure scheme; 


e include the influence of changes in the symmetric wind of the TC; 


e include the first-order effects of large-scale relative vorticity gradients of arbitrary 
magnitude and direction; and 


e include the first-order effects of moderate values of uniform relative vorticity of the 
large-scale environment. 


The close agreement between the TC propagation velocities predicted by this model 
and the equivalent numerical solutions verifics that the piecewise-analvtic construction 
technique is dynamically sound. The propagation vector errors from this model have 
small, well-defined biases that depend on the two adjustable parameters of the model in 
a predictable and dynamically reasonable manner. This property, combined with the 
capability of the model to predict TC propagation and gyre structure for many realistic 
combinations of TC structure and environmental windfields, suggests that this model has 
significant potential for use as either: i) an initialization scheme for barotropic numerical 
forecast models such as SANBAR; il) or as part of a “propagation + advection” track 
prediction aid. A demonstration of the potential usefulness of this model as an initial- 
ization tool is provided in Chapter V. 

The success of this theoretical mode! in reproducing the numerical results of CW 
and FE provides strong evidence that nonlinear self-advection, rather than hnear prop- 
agation. is the barotropic mechanism that contributes to the observed propagation of 
TC's. Itis important to emphasize that the individual effects of 6, environmental shear 


and environmental vorticity gradients are integrated within this self-advection model, 
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since those processes together determine the phase and amplitude of a wavenumber one 
interaction flowfield that advects the TC relative to the large-scale environment. By 
providing a unified theory for barotropic TC propagation, this model represents an im- 
portant step toward the development of a general theory of TC motion. For the sake 
of brevity and to emphasize the unifving aspect of self-advection theory, this model will 
hereafter be referred to as the Barotropic Self-advection Model (BSAM). 
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V. BAROTROPIC SELF-ADVECTION MODEL APPLICATION 


In this chapter, the usefulness and versatility of the BSAM will be demonstrated in 
three important areas. First, the issue of how accurately the outer wind strength must 
be measured to adequately account for propagation of TC’s will be addressed by using 
tangential wind profiles based on composited data as input. An alternate approach of 
using the BSAM to predict an “effective” outer wind strength based on past TC propa- 
gation will also be outlined. Second, a quantitative assessment of the extent to which 
the asymmetric interaction flow will be accounted for in steering flow calculations will 
be made based on the wavenumber gyre structure predicted by the BSAM for the com- 
posite data. Suggestions on how to isolate the environment and interaction flow com- 
ponents will also be included. Finally, the feasibility of using the BSAM to initialize 


barotropic numerical forecast models will be demonstrated. 


A. BSAM PREDICTIONS USING COMPOSITE DATA 
1. Preliminary analvsis 

In Chapter IV, the BSAM propagation predictions were based on piecewise- 
analytic TC wind profiles that closely approximated the exponential wind profile (4.22 
used in recent numerical studies (CW and FE). A similar exponential profile was also 
used by DeMaria (1985; DM). Exponential TC wind profiles were used in the BSAM 
to verify the accuracy of the theoretical predictions relative to equivalent numerical sol- 
utions. In light of the importance of TC structure on barotropic propagation, it 1s im- 
portant to examine how well exponential profiles approximate TC winds, and whether 
some other functional form might be better for propagation-prediction purposes. 

As shown in Fig. 5.1a, FE selected a “basic vortex” profile to have a radius of 
15 ms winds at 300 km to agree with typical composite observations of TC surface 
tangential winds, such as those given by Merrill (1984) for large Atlantic hurricanes. 
For comparison, a “composite pressure-averaged typhoon” was calculated by taking a 
950-150 mb pressure-weighted average of the western North Pacific composite TC 
tangential wind data of Frank (1977; his Fig. 9). The basic vortex has greater intensity 
and less outer wind strength than typical typhoons or large hurricanes, as represented 
by either the surface or the pressure-averaged tangential winds (Fig. 5.1a). The com- 
posite pressure-averaged wind profile for NW Pacific typhoons 1s significantly weaker 


at all radii than the composite surface wind profile for large Atlantic hurricanes because 
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Fig. 5.1 (a) Radial profiles of Fiorino and Elsberry’s “basic vortex” (dashed), the com- 
posite surface tangential winds of large Atlantic hurricanes (solid squares; Merrill 1984) 
and the tangential winds of the composite pressure-averaged typhoon (open squares, 
Prank 1977). Ilere ’,-=35 msand R,.= 100 km. (b) Analytic approximations to the 
radial prolile of the composite pressure-averaged typhoon (open squares) as defined by 
(4.22) with b= 0.6 (dashed) and (5.1) (solid; parameters im mset). dlere 1’,,.— 20 11s. 
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the tangential wind components decrease with elevation in TC’s. Using surface wind 
profiles to initialize barotropic models may cause an overestimation of TC propagation 
associated with barotropic processes. The composite pressure-averaged typhoon wind 
profile will be used here. 

Both FE and DM adjusted the maximum wind scale V’,,. and the b parameter 
in (4.22) to better fit the larger and weaker average typhoon;hurricane. Although the 
exponential profile can be adjusted to approximate the wind profile of the composite 
pressure-averaged typhoon (Fig. 5.1b), the exponential profile does not represent the 
curvature (and thus the vorticity gradient) of the composite profile nearly as well as the 


analytic function 
vs(r) = Dare SD — 1) + D3 =a) (5.1) 


which is a modified Rankine profile with a three-term Taylor series. Thus, (5.1) will be 
used as an analytic representation for composite TC tangential windfields below. 
2. Propagation speed versus composite TC strength 

By adjusting the curvature of the piecewise-analytic wind profile via the param- 
eter X, in (4.6), the outer wind strength of the composite pressure-averaged typhoon 
may be underestimated, closely approximated and overestimated in the 300-800 km 
annulus (Figs. 5.2-5.4 respectively), even though the inner wind speeds remain essenually 
unchanged. In each case, the interface radius R, is located at 300 km as in FE so that 
the same strength change methodology is used. An additional parameter C; is present 
(see insets) because the functional form of ry, in the Dispersion Region has been modified 


to 
Op ey Ce ee Cy (5.2) 


to require that €f,/¢cr be continuous at the transition radius. Although such a change 
adds two additional terms to (4.17b), it has a negligible influence on the BSAM-predicted 
propagation, and serves only to avoid overly large jumps in the piecewise-analytic 
vorticity gradient at R;. The magnitude of the analvtic vorticity gradient at R; is about 
four times smaller for the composite profile (Fig. 5.2b) than for the exponential profiles 
used earlier (e.g., Fig. 4.2b). Also, the change of vorticity gradient sign occurs well into 
the Dispersion Region for the composite data duc to the significantly greater outer wind 
strength. The vorticity gradient of the piecewise-analvtic profile results in some under- 


estimation of the analytic vorticity gradient at R,. Tlowever, it 1s not certain whether 
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Tig. 5.2 (a) Radial wind profiles of the composite pressure-averaged typhoon (dashed) 
as approximated by (5.1) im Fig. 5.1b and a piecewise-analytic profile (parameters in in- 
sect) that overestimates outer wind strength. (b) Piecewise-analytic and analytic vorticity 
eradient profiles corresponding to the winds in (a). In both cases, }’,.=20 m/s and 
Rye = 100 km. 
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Vig. 5.3 As in Vig. 5.2, except for a piecewise-analytic wind profile that closely approx- 
imates the outer wind strength of the composite pressure-averaged typhoon in the 300 
to 800 km radius “critical annulus.” 
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Vig. 5.4 As in Fig. 5.2, except for a piecewise-analytic wind profile that underestimates 
the outer wind strength of the composite pressure-averaged typhoon. 
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this will introduce a speed bias into BSAM propagation predictions because the coarse 
radial resolution and limited accuracy of the composite data allows for some subjectivity 


in selecting the curvature and fit of (5.1). 


Table 5.1 TC propagation velocities (columns 5,6) predicted by the BSAM for 
piecewise-analvtic Wind profiles (see Figs. 5.2a-5.4a) that underestimate, approximate 
and overestimate the outer wind strength of the composite pressure-averaged typhoon 
(Cases 1-3 respectively). In each case, a quiescent environment has been assumed. The 
wind speed at the transition radius R; (column 3) and at 550 km (column 4) are also 
shown. 


¥se(550 kin) 
(m s) 





The propagation speeds and directions predicted by the BSAM (Table 5.1) in- 
crease With increased TC outer wind strength, which is in general agreement with the 
results of Chapter IV and previous studies. Quasi-linear relationships exist between the 
propagation speed (column 5) and the wind speed at the transition radius (column 3), 
and the wind speed at 580 km (column 4). The first relationship has more relevance to 
the dvnamics of the BSAM since r,.(R;) is a rough measure of the amount the linear 
asymmetric forcing generated in the Dispersion Region, which acts to strengthen the 
wavenumber gvres and thus increase the propagation speed. The second relationship 1s 
also dynamically relevant because 550 km would be the center of FE’s “critical annulus” 
between 300 and 800 km. The 1 ms propagation speed increase for a 2.4 m’s increase 
in TC tangential winds at 550 km represents a potential for about 85 km of forecast error 
in 24h. This result suggests that TC outer wind strength must be measured to within 


+1 m s in order for numerical models to be initialized with sufficient accuracy to avoid 


significant forecast errors due to misrepresentation of TC strength. Important related 
questions are: 1) do numerical model heating and dissipative parameterizations tend to 
maintain the tangential wind profile to this accuracy throughout the forecast integration; 
and i1) are changes in TC strength that occur during a forecast period accurately re- 
produced by numerical forecast models? 

This sensitivity of TC propagation on outer wind strength poses a significant 
challenge for present observing system technologies. However, the method just used to 
demonstrate the wind-strength propagation-speed relationship offers a potential alter- 
nate approach to determine an “effective” TC strength. In Table 5.1, the independent 
adjustment of a single parameter X, (R; varies as a result) causes related changes to TC 
outer wind strength and propagation speed. It is therefore possible, in principle, to seek 
a value of X, such that the propagation speed predicted by the BSAM matches the ob- 
served Value at the initial time. The associated wavenumber one gvre structure predicted 
by the BSAM could then be used to initialize a barotropic forecast model like SANBAR, 
or provide guidance in selecting an appropriate bogus vortex for baroclinic forecast 
models. 

Iinplicit in such a scheme to predict TC strength are the requirements that: 
1) the large-scale environment is known sufficiently well to provide f,,., @ and Z;. for the 
BSAMI: 11) the track of the TC is resolved well enough to compute propagation relative 
to the barotropic steering of the environment; and it) an estimate of TC speed at about 
300 km is available. Since it is likely that many combinations of inner intensity and 
outer wind strength can give similar propagation velocities, such a scheme might not 
need to reproduce the outer wind structure of the TC, but only an inner- 
intensity outer-wind-strength combination (ie., an “effective” TC strength) that will 
produce the correct propagation. In addition, possible baroclinic contributions to TC 


propagation might also be parameterized to some extent by such an approach. 


B. INTERPRETATION OF COMPOSITE TC PROPAGATION VECTORS 

In Chapter I. comparison of observations of TC propagation with numerically de- 
rived depictions of f-induced wavenumber one gyre structure (FE) led to the preliminarv 
conclusion that the uniform flow region of the gyres would not be accounted for in 
steering flows calculated over an annulus of 5-7° lat. radius. However, the wavenumber 
one gvres for the composite pressure-averaged typhoon due to f-forcing only (Fig. 5.5) 
are somewhat larger than the gvres associated with the basic vortex profile of FE (Ing. 


4.4). The well resolved gyre structure provided by the BSAM will be used to calculate 
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Fig. 5.5 Wavenumber one gtre streamfunction (dashed, negative) associated with Table 
5.1 Case 2. Contour interval ts 2.0 x 10°171?/s. Dotted circles have radii of 2°, 4°, 6°, 8° 
and 10° lat. The middle circle corresponds to the center of the 5°-7° lat. radius annulus 
typically used to compute steering from composite data (e.g., Chan and Gray 1982; 
Holland 1984). The parameter R,,.= 100 km. 
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the extent to which TC self-advection is included in such steering flow calculations or is 
manifested as propagation. 

By averaging the velocities around various circles centered on the BSAM-predicted 
gyres (Fig. 3.5, dotted lines). a set of velocity vectors representing interaction flow con- 
tributions to environmental steering may be obtained (Fig. 5.6a). Subtracting such 
vectors from the BSAM-predicted propagation velocity for the composite pressure- 
averaged typhoon (Table 5.1: columns 5.6) gives a set of vectors that represent the 
propagation of the storm relative to the steering computed at a particular radius (Fig. 
5.6b). For composite steering computed at a 2° lat. radius, most of the interaction flow 
is accounted for in the steering, whereas the propagation vector is about 0.4 m s to the 
west. In contrast, at the 6° lat. radius where steering is typically computed, all of the 
westward component and about one quarter of the northward component are manifested 
as propagation. This result is consistent with the structure of the wavenumber one gvre 
uniform flow region in Fig. 5.5, which is oriented essentially north-south at 6° lat. radius, 
but still has a velocity nearly as large as the innermost portion of the flow. This result 
may then explain why TC propagation computed from composite data over a 5-7° 
annulus tends to have a predominantly westward orientation (Fig. 1.1). Although the 
analwsis here has used gvres associated with only the influence of planetary vorticity, the 
results would also apply to the more realistic environmental situations of Chapter IV 
Section D.3, if appropriate phase shifting of the vector patterns is taken into account. 

As discussed in Chapter III, the angular windshear in the symmetric TC produces 
a down-shear tilting of the inner portion of the wavenumber one gyre that barotropically 
stabilizes the TC to asvmmetric forcing. It is this down-shear tilt that produces the 
characteristic cyclonic rotation and increasing magnitude of TC propagation relative to 
steering flows computed at increasingly larger radi (Fig. 5.6b). This propagation vector 
pattern was evident in early compositing studies (George and Grav 1976: their Figs. 12 
and 13), and is a persistent feature of the recent compositing results based on much 
larger data sets presented bv W. Grav at a recent Office of Naval Research Workshop 
on Tropical Cyclone Motion (Ulsberrv 1989). Such vector patterns mav be interpreted 
as evidence that observed TC propagation is a manifestation of self-advection by 
wavenumber one gvres produced by asymmetric forcing. Although baroclinic processes 
mav also contribute to the total wavenumber one gvre pattern, the general agreement 
of observed propagation (Figs. 1.1) with the predictions of barotropic theory (Table 4.4) 


sugzesis that barotropic processes are the dominant influence. 
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Fig. 5.6 (a) BSAM-predicted propagation (C) for the composite pressure-averaged 
typhoon and contributions to environmental steering by the associated BSAM-generated 
wavenuinber one gyre flow evaluated around the 2,4,6,8 and 10° lat. radius circles of Fig. 
5.5. (b) Diflerences between composite pressure-averaged typhoon propagation (C) and 
the wavenumber one gyre contributions to steering of (a). 
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For the BSAM to be used as a numerical model initialization tool or as part of a 
propagation + advection forecasting aid, the interaction flow must be effectively excluded 
from the environmental windfield. Although filtering is often used in an attempt to re- 
move the influence of the 1C from the large-scale environment, alternate methods to 
compute steering that better exclude the uniform flow region of the gvyres should also 
be investigated. For example. the steering might be computed along circles as in Fig. 
5.5, but including only the 90° arcs centered on a line normal to the direction of the 
large-scale absolute vorticity gradient (Fig. 5.7). Using this approach (Fig. 5.8a), the 
interaction flow contribution to the 6° lat. radius steering accounts for only half the 
northward component of propagation, and the contribution to the 8° lat. radius steering 
is essentially zero. As a result, nearly all of the interaction flow influence is manifested 
as propagauon relative to steering computed at 8° lat. radius or larger (Fig. 5.8b). Thus, 
it mav be desirable to depart from traditional steering computation methods to better 


separate the large-scale environment from the influence of the TC. 


C. NUMERICAL MODEL INITIALIZATION 

The wavenumber one gvre structures provided bv the BSAM (based on cither an 
external specification or an internal prediction of TC propagation) should be ideal for 
initializing a nondivergent, barotropic (NDBT) numerical model so that quasi-steady 
motion occurs immediately. Comparisons with the NDBT numerical model of CW will 
be used to provide a preliminary verification of this assertion. The model uses a 
f-plane approximation, a 4000x4000 km channel domain with an east-west cyclic 
boundary condition and a resolution of 20 km. 

1. Quiescent environment predictions 

Using the symmetric, exponential TC wind profile of Table 4.2 Case 2, the CW 

model predicts a nearly constant f-induced propagation direction of about 330° during 
the first 48 h (Fig. 5.9). In contrast, the propagation speed increases from Zero initially 
to about 2.65 ms at 48 h. During this period, the asymmetric gyres develop. and then 
the propagation speed becomes almost constant (2.8 m's) after the gyres are established. 
In a real storm. these gvres are presumably present continually, and change only in re- 
sponse to Variations in TC structure or environnmental forcing. The “forecast error” 
associated with the transient acceleration period evident in Fig. 5.9 mav be estimated 
by comparing the 0-48 h track with a 48 h displacement between say 36 h and 84 h 


during which the propagation speed is quasi-steady. The comparison is facilitated by 
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Fig. 5.7 As in Fig. 5.5, except for including in steering flow computation only the winds 
in the 90° arcs to the right and left of the planetary vorticity gradient. 
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Fig. $.8 As in Fig. 5.6, except for vectors calculated based on the arcs shown in Fig. 5.7. 
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Fig. 5.9 f-induced propagation tracks and 6-hour positions predicted by the numerical 

model of Chan and Williams (1987) for an initial TC wind profile corresponding to Table 


4.2 Case 2 from 0-48 h (solid circles), and for the period 36-84 h (crossed circles) except 
translated so that the 36 h position corresponds to the initial position of the TC. 
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translating the 36-84 h track so that the 36 h position of the TC corresponds to the in- 
itial position. 

The along-track error of about 150 km suggested by this approach may be an 
underestimation since the TC has slowed to about 2.6 m's by 84 h due to a weakening 
of outer wind strength as the TC is displaced northward on a f-plane (cf., FE). The 
much smaller cross-track error is caused by the slightly more northward track of the TC 
from 36-84 h as compared to 0-48 h. Since this error is small, and may be associated 
with different TC strength during the two 48 h periods or boundary influences, it will 
be ignored. 

The numerical model is now modified to include the wavenumber one gyre 
Structure generated by the BSAM given the propagation velocity of 2.8 m’s at 330° 
(x= 120°) reported by CW, or the BSAM-predicted propagation velocity (Table 4.2, 
Case 2), as shown in Figs. 5.10a and 5.10b respectively. In both cases, essentially steady 
propagation of the model TC occurs immediately (Fig. 5.11). The propagation tracks 
for the two asymmetric initializations agree remarkably well, despite the 12° difference 
between the initial direction of propagation predicted bv the BSAM and the value ob- 
served by CW. This indicates that the numerical model has rapidly adjusted for differ- 
ences between the analvtically-predicted gvres and those that would be generated from 
an initially svmmetric TC wind profile (e.g., Fig. 4.4.). This is consistent with the results 
of Chapter III regarding short adjustment times near the radius of maximum winds, 
since only the’ central uniform flow region of the BSAM-generated gvres is altered by 
modest changes in propagation direction. Thus, the more important aspect of the in- 
itialization with the BSAM-generated gvres included is that the proper propagation 
speed is established immediately. 

Using the 36-84 h track in Fig. 5.11 as a benchmark, the two initializations of 
the numerical model with f-induced gyre structures produce an along-track forecast er- 
ror of about 60 km. As indicated above, this error is likely an overestimation, because 
the 36-84 h track includes an anomalously slow portion around $4 h. Thus, the along- 
track error at 48 h has been reduced from greater than 150 km to a value under 60 km. 
Part of the apparent cross-track error of 38 km relative to the translated 84 h position 
of the TC is an artifact of the more northward track of the TC from 36 h to 84h, as 
mentioned above. Comparing the 0-48 h track in Fig. 5.9 with the 0-48 h tracks in Tig. 


3.11 muicates that the cross-track error is actually about 20 kin. 
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Fig. 5.10 Wavenumber one gyre streamfunction patterns generated by the BSAM for 
the symmetric 1C of Table 4.2 2 Case 2 for (a) Chan and Williams propagation velocity 
of 2.8 m’s at 330° (#=120°). and (b) a BSAM-predicted propagation velocity of 2. 65 
ms at ¢=132°. Areal extent of the figures corresponds to the domain of the numerical 
model, and the streamfuncttons have been linearly adjusted to zero within 20 gridpoints 


of the domain boundaries. 





300 


(KM) 


1 
4 


Y~-DISTANCE 


200 





-500 -400 -300 200 -100 0 100 
X—DISTANCE (KM) 


Fig. §.1]| The f-induced propagation tracks and 6-hour positions predicted by the nu- 
merical model initialized with the symmetric TC as in Fig. 5.9, except including the gvre 
structures of Fig. 5.J!a (solid circles) and Fig. 5.106 (open circles). For comparison, the 
translated 36-§='h track (crossed circles) from Fig. 3.9 1s also shown. 


2. Results in a linearly-sheared environmental current 
As shown in Table 4.3, the BSAM predicts significant changes in TC propa- 
gation speed and direction for different values of uniform environmental vorticity. In- 
itialization tests may also be made with the CW numerical model by including a 


linearly-sheared zonal current defined bv 
Ups a Lee ys 3 (5.3) 


which has zero velocity at the initial position of the TC. The numerical model is then 
integrated for 48 h using the symmetric TC vortex of Table 4.2 Case 2 with and without 
the BSAM-predicted gvres for Table 4.3 Cases | and 2. Only BSAM-predicted propa- 
gation velocities can be used here since numerical predictions of TC propagation relative 
to a sheared zonal current are not presently available, although work is proceeding in 
this area (Willams 1989; personal communication). 

Comparison of the 0-48 h tracks of the initially symmetric and BSAM-initialized 
TC’s in the uniformly anticvclonic (Tig. 5.12) and cyclonic (Fig. 5.13) zonal currents 
reveals a peculiar aspect that warrants detailed analysis. In the anticvclonic shear case, 
the distance between the 24 h and 30h positions of the BSAM-initialized TC is equiv- 
alent to a speed of 3.3 ms. Since the TC is moving essentially due north during this 
period, such a movement can only be due to a f-induced propagation that is substan- 
tiallv faster than the 2.3 m s initially provided by the BSAM-generated gyre. This result 
may be due to the BSAM having a slow bias in the anticyclonic shear case. However, 
other evidence suggests the propagation speed increase may be associated with a non- 
steady process present in the numerical model. For example, the 42 h to 48 h movement 
of the BSAM-initialized TC is equivalent to 3.9 m’s, which is significantly faster than the 
24h to 30 h speed computed above. In additional, the 42-48 h speed of the initially 
symmetric TC in Fig. 5.12 is faster (3.0 ms) than the counterpart in a quiescent envi- 
ronment (Fig. 5.9), which was noted earlier to have a 42-48 h speed of 2.65 m/s. Thus, 
the evidence strongly suggests that the influence of uniform anticyclonic environmental 
shear causes a nonsteady alteration of f-induced TC propagation that cannot be ac- 
counted for by a steady-state model such as the BSAM. 

A similar result is not found for the cyclonic shear case. The 42-48 h speed of 
4.2 ms for the BSAM-initialized TC in cyclonic shear (Fig. 5.13) also is substantially 
faster than the propagation of 3.0 ms generated by the BSAM gvres (Table 4.3; Case 


2, column 3). However, since this track is largely westward, a significant fraction of this 
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Fig. 5.12 Tracks and 6-hour positions predicted by the numerical model for B-induced 
TC propagation in an zonal current with linear anticyclonic shear of 
Z,- = —5.0x 10s"! using the symmetric initial TC wind profile of Fig. 5.9 with (open 
circles) and without (solid circles) BSAM-generated wavenumber one gyres. 
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5.12, except for cyclonic shear of Z-- = 5.0% 10.*5* 
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total motion is due to the 1.5 m/s westward speed of the environmental current at that 
location. Subtracting the environmental velocity vector results in a B-induced propa- 
gation velocity of 2.9 m sina direction ox144°. The difference between this result and 
initial propagation specified in Table 4.3 Case 2 is insignificant compared to the large 
change in propagation velocity that occurs in the anticyclonic case above. 

The solutions of Chan and Williams (1989) for B-induced TC propagation in a 
linearlv-sheared zonal current also contain anomalously fast propagation in anticyclonic 
shear, but no measurable changes in cyclonic shear. Thus, this unexplained behavior is 
not limited only to the BSAM. The increased propagation speed in anticyclonic shear 
is certainly not intuitively plausible, since the results of FE (if applicable) would suggest 
that the weakening of the outer wind strength of the TC due to a superposition of an 
anticyclonically sheared current should cause slower propagation. Additional numerical 
integrations are being performed to gain a better understanding of the dynamical or 
perhaps numerical basis for this behavior (personal communication, R. T. Williams 
1989). 

The potentially nonsteady influence of linear environmental windshear on 
B-induced TC propagation make a precise determination of how well the BSAM has in- 
itialized the CW numerical model difficult. However, several qualitative observations 
and comments can be inade. It is evident in Figs. 5.12 and 5.13 that significant cross- 
track as well as along-track errors exist between the positions of initially symmetric and 
BSAM-initialized TC’s when a spatially variable environment is present. As in the 
quiescent environment case above, the along-track differences are due to the nonzero 
initial speed of the BSAM-initialized TC. The cross-track differences are in part due to 
a slightly different starting direction for the BSAM-initialized TC, which in effect ac- 
counts for phase shifting of the wavenumber one gvres by the “past” influence of envi- 
ronmental shear. Additional research will be necessary to understand to what extent the 
present asymmetric structure of TC’s reflects past environmental influences. The cross- 
track differences are also due to the nonzero initial propagation speed of the 
BSAM-initialized TC which. continually subjects the TC to different environmental 
winds compared to the initially symmetric TC. Thus, the tracks of the BSAM-initialized 
and initially symmetric TC diverge with time. The divergence is much greater in the 
cvclonic environment case due to the greater speed of propagation and the cyclonic ro- 
tation of the f-induced gvres such that a larger zonal motion component (propagation 
~ environmental advection) is produced. It should also be emphasized that the 48 h 


along-track and cross-track differences of up to 280 km and 80 km respectively in Figs. 


to 
~] 


5.12 and 5.13 indicate that initializing a numerical model with only a symmetric TC wind 
profile has much more serious consequences when spatially spatially variable environ- 


mental winds are present than for a quiescent environment (Fig. 5.11). 


D. SUMMARY 

In this chapter, three basic but important applications of the BSAM have been 
demonstrated. First, the BSAM has been used with composite observations to predict 
the barotropic propagation due to f for a composite pressure-averaged typhoon, and to 
predict the change in propagation speed that nught be expected for a typhoon with 
larger or smaller outer wind strength. The results suggest that misrepresenting the outer 
wind strength by about 2.5 m s will cause a | m’s error in barotropic propagation, which 
represents a potential forecast error of 85 km in 24h. This application also outlined a 
method by which the BSAM might be emploved to predict an “effective” outer wind 
strength based on a priori knowledge of TC propagation. 

Second, analysis of the f-induced gyre structure predicted by the BSAM for the 
composite pressure-averaged typhoon provided evidence that the TC propagation vec- 
tors in composite data (Fig. 1.1) are indeed manifestations of wavenumber one gyres 
associated with self-advection processes. An alternate method for computing environ- 
mental steering that better separates the large-scale environment from the self-advection 
flow was also demonstrated. 

Finally, and most importantly, the potential usefulness of the BSAM as an initial- 
ization tool for barotropic numerical forecast models was demonstrated. In the case of 
a quiescent environment. initializing the numerical model with BSAM-predicted 
wavenumber one gvres reduced the along-track forecast error at 48 h from more than 
150 km to less than 60 km. When a linearlv-sheared zonal environmental current was 
included, significantly larger along-track and cross-track differences of up to 280 km and 
SO km respectively were observed between the 48 h positions of an initially symmetric 
TC and a BSAM-initialized TC. 
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VI. CONCLUSION 


The principal achievement of this research has been the development, testing and 
preliminary application of a comprehensive theoretical model for predicting tropical 
cyclone (TC) propagation and associated asymmetries based on the barotropic self- 
advection process identified in recent numerical studies. The model is comprehensive in 
the sense that it includes the first-order effects of all of the dynamical processes that are 
presently understood to be important to barotropic TC propagation: gradients of plan- 
etary and environmental vorticity. changes in TC wind structure, and environmental 
windshear. Since the model is based on a single self-advection principle through which 
individual external forcing effects collectively determine the phase and amplitude of a 
Wavenumber one gvre flow that advects the TC, the nomenclature Barotropic Self- 
Advection Model (BSAM) has been adopted. Because of the dynamical complexity of 
barotropic TC self-advection, a rather lengthy chain of hypotheses, assumptions and 
approximations were emploved to make the problem analsticallv tractable. The princi- 
pal ‘links’ in this development chain are summarized in Section A and a summarv of 


results follows in Section B. 


A. OVERVIEW OF PRINCIPAL ANALYSIS TECHNIQUES 
1. Dissection of the problem and model formulation 

Based: on prior numerical findings that TC propagation relative to steering re- 
sults from wavenumber one asymmetries induced by TC-environment interaction, the 
total flow field has been partitioned into three components: 1) a specified large-scale en- 
vironment: il) a s¥mmetric cyclone circulation; ili) and an unknown asvmmetic inter- 
action flow component that is presumed to be responsible for TC propagation. A 
transformation to a reference frame moving with the TC has been utilized, and the 
symmetric TC component is assumed steady in such a reference frame to further simplif¥ 
the problem. By subtracting from the total problem an equation governing the large- 
scale environment in the absence of the TC, an equation for the evolution of the inter- 
action flow in response to external forcing has been obtained. Thus, the partinoning and 
transformation process significantly simplified the problem by focusing the analysis on 
that portion of total flow field that is relevant to TC propagation. 

The mathematical complexity associated with the polar symmetry of the prob- 


Jem and the extreme radia! anisotropy of the symmetnc TC flow vaniables makes it very 


difficult to find a single solution that is analytic in the whole model domain. Thus, it 
has been assumed that a piecewise-analvtic solution obtained using a matched- 
asymptotic approach will be a reasonable approximation. In the preliminary model de- 
velopment (Chapter II), this assumption entailed dividing the model domain relative to 
a “transition radius.” Inside this radius is a Self-advection Region in which mutual 
advections by the symmetric TC and interaction flow dominate, and outside this radius 
is a Dispersion Region in which those advections are unimportant relative to the ex- 
ternal asymmetric forcing. After formally showing that only wavenumber one processes 
can contribute significantly to TC propagation, the problem has been further simplified 
by retaining only the first-order contributions to wavenumber one processes in the 
interaction flow governing equation. 

To overcome the problem of radial anisotropy noted above, the model domain 
has been subdivided further (Chapter IV) by assuming that the symmetric TC winds in 
the Self-advection Region can be reasonably represented by a set of piecewise-analytic, 
modified Rankine segments. A solution valid for the Whole model domain then has been 
assembled from solutions valid in individual annular regions by imposing matching 
conditions on solutions at the annulae interfaces. Thus, the intractability of the full 
problem has been overcome bv subdividing the problem sufficiently to make the indi- 
vidual pieces analytically tractable. The success of this systematic dissection approach 
to modeling TC self-advection has been demonstrated by several comparisons with 
complete numerical solutions and with some observational results. 

2. Free versus forced transient analysis 

An examination of the response of initially imposed perturbations on a Rankine 
vortex provides important theoretical insights into why TC-like vortices in the 
barotropic numerical models propagate steadilv in a slightly deformed state in response 
to persistent environmental forcing, rather than rapidly distorting and dispersing. The 
unforced initial value problem analyzed in Chapter III represents the homogeneous 
counterpart to the forced steady propagation problem addressed in Chapter IV. By 
choosing initial conditions with the same spatial structure as in the forcing terms of the 
steady-state problem, the forced transient adjustments in numerical TC models may in- 
terpreted in terms of the analytically tractable unforced problem. This modeling tech- 
nique permits a formal mathematical analysis of barotropic vortex adjustment to 


unposed asvmunetries that was not possible for the related forced problem. 
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3. Closing the Barotropic Self-advection Model (BSAM) 

Although the transformation to a reference frame moving with the TC was an 
important step in dividing the TC self-advection problem into manageable pieces, the 
presence of the propagation speed and direction in the interaction flow equation re- 
presented two additional unknowns. Thus, a closure scheme based on numerical results 
has been devised in which it is hypothesized that the correct propagation velocity will 
be associated with a particular wavenumber one gyre structure that has an approximate 
balance between vorticity tendencies due to the advection of large-scale absolute 
vorucity by the symmetric TC and the shearing of interaction flow vorticity by the 
symmetric TC winds. Specifically, the advection of the vorticity of the symmetric TC 
by the wavenumber one interaction flow is assumed to account for all of the propagation 
vector, except for a 0.4 m s component associated with linear Rossby wave propagation 
of the dominant scale of the TC. The closure scheme depends sensitively, but predict- 
ably, on the two adjustable features of the BSAM associated with the transition radius. 
In a comparison with equivalent numerical model solutions, the closure scheme error 


was 0.25 ms and 12° without any adjustment of the two available parameters. 


B. SUMMARY OF RESULTS 
1. Barotropic vortex stability 

The first important result of this research is that the angular windshear of the 
symmetric circulation of TC-like vortices acts to make the vortex barotropically stable 
to small asymmetric perturbations from axisymmetry. The windshear in the symmetric 
circulation tilts the perturbations downshear. which results in a barotropically stable 
transfer of kinetic energy from the perturbation to the “basic state” represented by the 
symmetric vortex. The damping dependence is an algebraic continuous spectrun)-type 
response analogous to the damping of perturbations imposed on an f-plane Couette 
flow. The damping rate is proportional to perturbation wavenumber, and is faster 
(slower) for perturbations with a radial structure that decays rapidly (slowly) with in- 
creasing radius. A wavenumber one perturbation associated with the advection of 
planetary vorticity by the symmetric TC loses about 80% of the imtial kinetic energy 
by 24 h. which agrees well with the adjustment times in recent numerical modeling 
studies of TC adjustment on a f-plane. Thus, the quasi-steady propagation of TC’s in 
barotropic models may be regarded as a balance between the transfer of kinetic energy 


fron) the symmetric TC to the wavenumber one interaction flow due to external asymi- 
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metric forcing and the transfer of kinetic energy to the symmetric TC from the 
wavenumber one asvmmetry due to the symmetric angular windshear. 

The presence of significant tangential windshear in TC’s, and the qualitative 
similarity of TC propagation tracks in baroclinic and barotropic models, strongly sug- 
gests that this barotropic vortex stability mechanism also may play a fundamental role 
in actual TC propagation. It should be emphasized that the asymmetry-damping influ- 
ence of symmetric angular windshear is not limited to vortices scaled to approximate 
TC’s. Rather, the shear-damping mechanism is applicable to any vortical flow that may 
be reasonably approximated by barotropic dynamics, and in which the strength of the 
vortical flow is strong compared to other influences over a significant horizontal area. 
Examples are Gulf Stream rings and other intense ocean eddies. 

2. Dependence of propagation speed on TC strength 

The piecewise-analytic symmetric TC wind profile in the BSAM permits propa- 
gation predictions for various tangential wind profiles. During testing of the closure 
scheme, it has been demonstrated that the BSAM-predicted propagation speed depends 
strongly on the magnitude of the symmetric vorticity gradient in the vicinity of the 
transition radius. This result complements the conclusions of Fiorino and Elsberrv 
(1989) regarding strength changes quite well, since the transition radius is typically lo- 
cated near the middle of their 300-800 km critical annulus. Specifically, the BSAM 
demonstrates that the speed of TC propagation depends on the second derivative of the 
TC tangential winds, which explains why the propagation; wind-strength dependence is 
so strong. 

The model has also been used with composite TC tangential wind profiles to 
obtain a practical estimate of the wind measuring accuracy required to avoid significant 
forecast error due to misrepresenting TC strength. Piecewise-analvtic wind profiles have 
been constructed that underestimate, approximate, or overestimate the outer wind 
strength of a composite pressure-averaged tvphoon. A quasi-linear relationship is found 
between the symmetric TC windspeed at 550 km and the associated propagation speed. 
Specifically, the speed of propagation increases by | m’s for a 2.4 m’s increase in v,.(550 
km), which suggests that the TC outer wind strength must be measured to about +1 
ms to avoid having 24 h forecast errors greater than 85 km. The capability to alter 
outer wind strength by changing only one parameter suggests that the BSAM might be 
used to estimate an “effective” TC outer wind strength given sufficiently accurate infor- 
mation about the environmental windfield, the past propagation of the TC relative to 


environmental steering and the intensity of the storm. 


3. Dependence of propagation on environmental vorticity 

The influence of uniform environmental vorticity in altering the phase of the 
wavenumber one gvres in the Dispersion Region can be included in the BSAM for values 
of Z,. (x 10-§ s-') inside the range -0.5 to 1.5. An increased (decreased) propagation 
speed and a more (less) westward propagation direction is predicted when cyclonic 
(anticyclonic) uniform vorticity is present compared to f-induced propagation in a 
quiescent environment. Although the direction changes are generally consistent with 
propagation vectors in an unpublished numerical study by Chan and Williams (1989), a 
verification for the BSAM-predicted speed changes is not available. An unresolved issue 
here is that a steadily increasingly propagation speed is evident in the Chan and Williams 
solutions for anticyclonic shear, whereas an essentially steady propagation speed is 
found for cyclonic shear. If the steadily evolving propagation speed and asymmetry 
situation is relevant to TC’s, then the question of how long the TC has been in a par- 
ticular environment may be important to track forecasting. 

BSAM predictions of TC propagation have been made for situations involving 
zonal and meridional gradients of environmental vorticity, the influence of 8, and a 
uniform component of environmental vorticity. For example, a TC north (south) of the 
subtropical ridge where a westward gradient of environmental vorticity is also present is 
predicted to have a more northward (westward) and faster (slower) propagation in gen- 
eral agreement with composite observations in the western North Pacific (e.g., Fig. t.1). 

4. Numerical model initialization 

An important potential application of the BSAM is to provide the initially 
symmetric bogus vortex of a barotropic numerical model with the correct wavenumber 
gvre structure so that quasi-steady propagation occurs immediately, and thus overcome 
a Well-known slow bias. Such an application may also be regarded as an indirect test 
of the validity of the modeling concepts, techniques and assumptions employed in de- 
veloping the BSAM. For a quiescent environment, the BSAM-predicted wavenumber 
one gyre produces virtually steady propagation regardless of whether an internally or 
externally-derived propagation velocity is used. Including the BSAM-generated gvres in 
the numerical model initialization reduces the along-track forecast error from more than 
150 km to less than 60 km. This preliminary result shows that the BSAM has significant 
potential as an initialization tool for operational barotropic forecast models such as 
SANBAR. 

The unresolved issue regarding nonsteady f-induced propagation in anticyclonic 


sheur mentioned above precludes a precise estimation of the potential 48 h forecast error 
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reduction that might be achieved by a numerical model initialization with the BSAM. 
Nevertheless, the along-track (cross-track) differences of up to 280 km (80 km) between 
the 48 h positions of TC with and without BSAM-initialization illustrates the potential 
sensitivity of TC track prediction to the initial asymmetric structure of the TC when a 
spatially variable environment is present. Using the adjustable parameters in the BSAM 
to minimize barotropic forecast model track error in some statistical sense may be a Vi- 
able approach for developing an initialization scheme that can improve barotropic nu- 
merical model forecast skill, and provide additional insights into the interaction of TC’s 


with the surrounding environment. 


APPENDIX A. COMPOSITE DATA CONVERSION 


Composite studies have generally characterized TC propagation in relative terms 
(e.g.. speed and direction differences) using a rotated coordinate system aligned with 
storm motion. This compositing methodology tends to make theoretical interpretation 
of the data difficult because the analytical and numerical studies predict that TC prop- 
agation will possess a particular orientation with respect to the direction of the large- 
scale vorticity gradient. In particular, a rotated storm-relative coordinate system would 
tend to obscure TC propagation associated with the gradient of the Coriolis parameter, 
since that gradient has a storm and environment-independent northward orientation. 
Thus, part of the difficulty in comparing theory with the composite observations may 
be readily overcome by representing TC propagation as a vector quantity in a north- 
oriented earth-relative coordinate system. 

The TC motion and steering flow composite data are taken from the latitude, di- 
rection, speed and intensity stratifications of Chan and Gray (1982) and George and 
Gray (1976) for the western North Pacific region, and the direction and recurvature 
stratifications of Holland (1984) for the Australian-Southwest Pacific region. Holland 
used a single steering flow definition based on a 800 to 300 mb pressure-weighted mean 
wind averaged over an annulus extending 5° to 7° lat. from the TC center. Although 
Chan and Grav used the same horizontal domain, several vertical averaging schemes 
Were tested. Only the Chan and Gray steering flow based on a surface to 300 mb vertical 
average is used here, since it most closely approximates the steering flow definition used 
by Holland. Vertically-averaged steering flows have been chosen for this analysis rather 
than individual steering levels (e.g., George and Gray 1976) to more appropriately com- 
pare the observations with the theoretical modeling results that are predominantly based 
on barotropic dynamics. 

Table A.1 summarizes the conversion process for the western North Pacific com- 
posite data. The columns labelled V,y, DD. V_ and D, contain previously published 


data. and the last two columns have been computed using the relationships 


veo a 
pe Ay 
Vp ean it) 


Table A.1 Original and converted composite TC motion data for the western North 
Pacific region. Column heading meanings: Vpy, is the speed of the steering flow compo- 
nent parallel to the direction of the composite TC minus the speed of the TC; DD is the 
difference between the direction of TC motion and the steering flow; V. and D, are the 
specd and direction of motion of the TC respectively; and V, and D, are the speed and 
direction of the steering flow respectively. The data in columns V,,, and DD are taken 
directly from Chan and Grav (1982). and the data in columns V;. and D, are taken di- 
rectly from George and Grav (1976). The data in columns V, and D, have been com- 
puted as described in the text. Directions are measured clockwise from North and the 
data in the last four columns are relative to a reference frame fixed to the surface of the 
earth. 


Composite Maer Eos 
Stratification (m s) ees ne S) ae pie (deg) 


Latitude: 


S20 SGT) 

<20F \ (LE) 
Direction: 

Westward (W) 


Northward (N) 
Eastward (E) 


Speed: 


Slon Ss) 
Moderate (MM) 
Fastih) 


Intensity: 
Weak (WK) 


Intense (1) 
Very Intense (VI) 





All column labels are defined in the table caption. Table A.2 1s analogous to Table A.1 


for the Austrahan-Southwest Pacific composite data. The columns labelled SD, DD. V5, 
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and D, contain previously published data, and the columns labelled V. and D;, have 
been computed using (A.1) and 


Ve — Ve eas SD . (A.3) 


The TC propagation vectors shown in Fig. 1.1] were then computed for each composite 


stratification using the last four entries in each row of Tables A.1 and A.2. 


Table A.2. Analogous to Table A.1 for the Australian-Southwest Pacific region. The 
column headings DD, V,.. Dc. Vg and D, have the same meanings as in Table A.1 and 
SD 1s the speed difference between the composite TC and steering. The data in columns 
SD, DD, V, and D, are taken directly from Holland (1984), except that the steering flow 
directions are measured clockwise from North. The data in columns V,. and D; have 
been computed as described in the text. 


Composite 
Stratification 


Direction: 


Westward (W) 
Southwestward (SW) 


Southward (S) 
Southeastward (SE) 
Eastward (E) 


Recurvature: 
Before (B) 


Near () 
After (A) 





APPENDIX B. PIECEWISE-ANALYTIC VORTEX CONSTRUCTION 


The following outline describes the procedure for constructing piecewise-analytic 
radial profiles of TC tangential wind and relative vorticity. The approach here is to use 
the linear degrees of freedom in (4.6), (4.16) and (5.2) to develop closed-form expressions 
for the coefficients A,, A,, A;, B,, B,, B; and C, such that the piecewise-analvtic profiles 
must have continuous windspeed and relative vorticity. The remaining nonlinear degrees 
of freedom are then used to adjust the fit of the piecewise-analytic profiles to the corre- 
sponding analytic profiles by selecting various values for X, and X, and assessing the 
results via interactive computer graphics. 


1. Given the trial values of R, and R, from Step 1 of the model solution procedure 
(Chapter IV Section A.4). choose trial values for v,(R,) and ¥,(R,). Generally R, 
will be in the vicinity of the radius of maximum winds, and thus y,(R,)=1 will pro- 
vide a good first estimate. The value of y.(R,) for the piecewise-analytic vortex is 
chosen to closely approximate the windspeed of the analvtic profile at R,. 


to 


Determine the coefficients A, and B, for (4.6) in the first annulus, by solving the 
2x2 linear system that results from applying the boundary conditions of Step 1 
above to (4.6). 


Loe } 


Determine the coefficients A, and B, bv requiring the windspeed and vorticity of the 
piecewise-analvtic TC structure to be continuous at R,. As in Step 2 above, this 
entails solving a 2x2 linear system. 


4. Repeat the process in Step 3 to determine the coefficients A, and B, in (4.16) by 
requiring the piecewise-analyvtic windspeed and vorticity to match at the transition 
radius R;, For the piecewise-analvtic profiles used in Chapter V, this step is mod- 
ified to include evaluation of C, in (5.2) by solving the linear 3x3 system that results 
from requiring the piecewise-analytic windspeed, vorticity and vorticity gradient to 
be continuous at R,. 


5. Using interactive computer graphics, vary the parameters X, and X, to adjust the 
piecewise-analytic windspeed and vorticity profiles to approximate the analyuc 
counterparts as described in Chapter IV Section B.1. During this step, also adjust 
R,. Ry. vs(Ro) and v.(R,) if necessary to achieve an acceptable fit. 
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